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Preface 

The  purpose  of  this  research  is  to  develop  an  analytical  model  for 
scheduling  check-out  servers  at  the  commissary.  The  goal  of  the  model  is 
to  insure  that  the  average  customer  waiting  time  remains  constant 
throughout  the  scheduling  period.  Such  a  model  can  save  commissary 
management  much  time  and  effort,  and  it  can  save  money  by  improving 
the  utilization  of  the  commissary  workers.  The  model  was  developed  to  be 
general  enough  to  be  used  at  any  commissary  in  the  Air  Force,  and  it  can 
also  be  used  at  many  other  different  service  organizations. 

1  want  to  extend  my  sincere  thanks  to  my  thesis  advisor,  Major 
Joseph  Lltko,  for  all  of  the  help  he  gave  me  in  the  development  of  this 
thesis.  Thanks  also  to  my  thesis  reader,  Dr.  James  Chrissis,  who  also 
provided  help  along  the  way.  The  biggest  thanks,  however,  go  to  my  wife 
Aneita  and  daughter  Elizabeth,  for  being  understanding  when  I  did  not 
give  them  the  time  that  they  deserved  during  the  past  18  months. 
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\  i .  ,  Abstract 

The  purpose  of  this  reocoreh  was  to  develop  an  analytical  model  that 
would  optimally  schedule  commissary  checkers  so  that  the  expected 
customer  waiting-time  would  remain  relatively  constant  throughout  the 
scheduling  period.  A  two-phase  model  was  developed  to  solve  the 
problem.  The  first  phase  of  the  model  used  dynamic  programming  to  find 
the  optimal  number  of  checkers  required  throughout  each  day  to  meet  the 
desired  customer  waiting-time  goal.  Since  checkers  cannot  be  scheduled  to 
work  arbitrarily  short  tours  of  duty,  a  second  phase  was  needed  in  the 
model  to  find  the  optimal  number  of  checkers  to  assign  to  allowable  shifts 
m  order  to  meet  the  optimal  requirements  determined  in  phase  one. 

A  simulation  was  developed  to  validate  the  checker  scheduling 
model.  It  was  found  that  the  scheduling  model  produced  acceptable  results 
until  the  last  few  periods  of  the  day.  Additional  servers  needed  to  be 
added  heuristically  near  the  end  of  each  day  to  obtain  the  desired  customer 
waiting  times,  '"<’r  i  s  ;  C  btcV'.r  i/Arj  e-bjii-'h*' 

Several  extensions  of  this  work  are  possible.  First,  an  improved 
approximation  for  customer  line  lengths  could  be  used  at  the  end  of  each 
day.  Use  of  such  an  approximation  could  eliminate  the  need  for  heuristic 
rules  in  scheduling  servers  during  the  last  few  periods  of  each  day. 

Second,  the  scheduling  algorithm  that  was  developed  did  not  account  for 
checker  lunch  breaks.  Accounting  for  lunch  breaks  complicates  the 
problem,  but  two  different  approaches  were  suggested  for  a  solution 
allowing  for  checker  lunch  breaks.  Finally,  a  third  phase  could  be  added  to 
the  model  that  would  allow  assignment  of  actual  workers  to  the  optimal 
shifts  determined  in  the  second  phase. 
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OPTIMAL  SERVER  SCHEDULING  TO  MAINTAIN 
CONSTANT  CUSTOMER  WAITING  TIMES 


I,  introduction 


Background 

Controlling  customer  waiting  times  at  service  organizations  such  as 
the  commissary  is  a  difficult  task.  Long  lines  are  a  common  cause  of 
customer  complaints.  At  the  commissary,  long  lines  do  not  usually  cause 
customers  to  leave  the  store.  However,  they  can  result  in  loss  of  future 
commissary  sales  by  causing  customers  to  do  future  grocery  shopping  at 
off-base  establishments.  Long  lines  can  also  impair  the  efficiency  of 
commissary  service  by  creating  congestion  in  the  aisles. 

Long  lines  are  not  the  only  problem  facing  commissary  management. 
The  other  extreme,  lines  that  are  too  short,  is  also  a  problem.  From  a 
customer’s  standpoint,  short  lines  are  ideal,  but  short  lines  cost  the 
commissary  extra  money.  To  attain  short  lines,  the  commissary  must 
employ  extra  check-out  servers  (checkers).  Since  each  store  is  only 
allocated  a  set  number  of  checker- hours  each  month,  a  store  may  not  have 
the  needed  checker-hours  available  to  achieve  short  lines.  Somewhere, 
between  long  and  short  lines,  an  ideal  exists.  Keeping  the  line  length  and 
the  corresponding  customer  waiting  time  at  this  ideal  is  difficult,  especially 
in  the  face  of  limited  total  monthly  checker-hours. 
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The  easiest  and  most  common  way  to  control  a  queue's  length  is  by 
varying  the  number  of  servers.  Obviously,  with  more  servers,  shorter 
lines  would  be  expected.  Analytical  expressions  relating  the  number  of 
servers  to  the  number  of  customers  in  line  are  commonly  available  as  long 
as  certain  assumptions  are  met.  One  of  these  assumptions  is  that  the  mean 
customer  arrival  rate  remain  constant.  Unfortunately,  the  mean  customer 
arrival  rate  at  the  commissary  (and  at  many  other  service  organizations) 
varies  throughout  the  day,  making  the  scheduling  of  servers  more  difficult. 

Specific  Problem 

Current  scheduling  of  commissary  checkers  requires  a  considerable 
amount  of  the  store  management's  time.  At  the  larger  stores,  the  time 
spent  on  scheduling  is  estimated  between  eight  and  fourteen  hours  per 
week  (Polk,  1988).  Efficient  allocation  of  servers  depends  mainly  upon  the 
experience  of  the  scheduler.  A  reliable  and  automated  method  for 
scheduling  the  checkers  would  save  management  time  and  potentially  save 
money  through  improved  efficiency. 

Research  Objective 

The  primary  objective  of  this  research  is  to  develop  an  analytical 
model  that  will  optimally  schedule  commissary  checkers  so  that  the 
expected  customer  waiting  time  is  constant  thi  oughout  the  day.  Some 
suhobjectives  associated  with  the  primary  objective  are: 

1.  Make  the  model  sufficiently  general  so  that  it  can  be 
used  at  any  Air  Force  commissary. 

2.  Specify  a  goal  for  the  mean  customer  waiting  time  and 
achieve  that  goal  throughout  the  month. 
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3.  Keep  the  total  number  of  checker-hours  scheduled 
during  the  month  below  a  given  maximum  number  of 
allowable  checker- hours  for  the  month. 

4.  Validate  the  analytical  model  using  simulation. 

5.  Observe  any  trends  in  server  requirements  during  a 
month. 

S_CQpfi 

The  monthly  scheduling  of  checkers  is  a  three  phase  problem.  In 
Phase  I,  checker  requirements  must  be  determined  for  the  scheduling 
period.  That  is,  the  ideal  number  of  checkers  required  to  achieve  the  mean 
customer  waiting  time  goal  are  found  in  Phase  1.  Then,  in  Phase  II,  all 
possible  checker  shifts  are  enumerated,  and  the  optimal  number  of  checkers 
are  assigned  to  each  shift  to  meet  the  requirements  calculated  in  Phase  1. 
Finally,  in  Phase  III,  actual  checkers  are  matched  to  the  optimal  shifts. 

This  research  effort  concentrates  on  Phases  I  and  11.  Phase  III,  which  is 
essentially  an  assignment  problem,  is  left  for  future  work. 

Plan  of  the  Report 

Chapter  1  has  introduced  the  problem.  The  background,  specific 
problem,  research  objective,  and  scope  of  the  research  were  discussed.  In 
Chapter  II,  the  literature  pertaining  to  this  problem  is  reviewed.  Included 
in  the  literature  review  are  sections  discussing  dynamic  programming  with 
resource  constraints,  fluid  approximations  to  queues,  manpower  shift 
scheduling,  integer  and  network  programming,  and  application  of 
lagrangian  relaxation  to  integer  programming.  Chapter  III  documents  the 
development  of  the  checker  scheduling  algorithm.  In  Chapter  IV,  a 
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simulation  is  used  to  validate  the  checker  scheduling  algorithm.  Finally, 
the  results,  conclusions,  and  recommendations  are  made  in  Chapter  V. 


1L  Literature  Review 


Overview 

There  are  six  main  areas  addressed  in  the  literature  review.  The 
first  section  reviews  dynamic  programming  under  constraints.  The  second 
section  discusses  a  simple  queuing  approximation  used  to  estimate  customer 
line  lengths.  In  the  third  section,  various  shift  scheduling  methods  are 
discussed.  The  fourth  section  of  the  literature  review  shows  how  certain 
integer  programming  problems  can  be  transformed  into  network 
programming  problems  and  the  fifth  section  outlines  how  these  network 
problems  can  be  solved.  The  final  section  of  the  literature  review  discusses 
the  lagrangian  relaxation  technique  for  solving  integer  programming 
problems  that  have  a  special  structure. 

Constrained  Dynamic  Programming 

Dynamic  programming  is  defined  by  Hillier  and  Lieberman  as  “...  a 
useful  mathematical  technique  for  making  a  sequence  of  interrelated 
decisions.  It  provides  a  systematic  procedure  for  determining  the 
combination  of  decisions  that  maximizes  overall  effectiveness  (Hillier  and 
Lieberman,  1988:332).”  If  a  problem  can  be  easily  split  into  stages  then 
dynamic  programming  should  be  considered  as  a  possible  solution 
technique.  At  each  stage  the  system  can  be  in  one  of  a  number  of  different 
states.  A  decision  is  made  at  the  present  stage.  The  effect  of  this  decision 
is  to  transform  the  system  state  at  the  present  stage  into  a  system  state  at 
the  next  stage.  The  mechanics  of  this  transformation  are  usually  defined 
by  a  transition  equation.  A  recursive  function  f  is  used  so  that  the 
decisions  made  at  each  stage  are  optimal.  The  difficulty  in  applying 
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dynamic  programming  is  in  defining  the  recursive  function  f  and  the 
transition  equation,  which  together  define  how  to  move  from  one  stage  to 
the  next.  Figure  1  shows  a  graphical  representation  of  dynamic 
programming. 
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Figure  1  Graphical  Representation  of  Dynamic  Programming 


One  typical  definition  of  the  recursive  function  f  at  stage  n  is 
(Denardo,  1982:162): 


max^min)  £  R*(n) + fn+i  (j)  J  for  n  <  N 
max^min)  £  R  *  (N)  J  for  n*  N 


where 


n  =  current  stage 
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i  *  system  state  at  current  stage 
j  *  system  state  at  next  stage 
N  »  total  number  of  stages  in  the  problem 
k  *  decision  made  at  stage  n 

R  *(n)  »  some  return  function  given  decision  k,  state  i,  and  stage  n 

Typically  the  return  function  is  some  cost  function  and  the  objective  is 
minimization.  For  this  case  R*(n)  would  represent  the  cost  at  stage  n 

of  decision  k ,  and  fn+l(j)  would  represent  the  total  cost  of  the  remaining 
stages.  Thus  the  current  decision  directly  affects  costs  at  the  current  stage 
and  indirectly  affects  costs  later  by  determining  the  next  system  state. 

A  dynamic  programming  problem  can  be  solved  by  starting  at  the 
final  stage  N  and  working  backward.  The  recursive  function  f  is 
calculated  for  every  possible  state  at  stage  N  .  To  calculate  f  for  a  given 
state  i ,  R^(N)  must  be  calculated  for  every  possible  decision  k  .  The 
recursive  function  f(i)  is  then  equal  to  the  minimum  value  of  Rj(N)  if  the 
objective  is  minimization.  Once  has  been  calculated  for  all  possible 
states  of  stage  N  ,  a  similar  process  is  followed  for  stage  N-l.  However, 
now  the  recursion  function  is  calculated  using  the  first  part  of  Eq  (1): 

rN_l(l)»  min  [R^(N-l)ff«0)]  (2 

One  difficulty  in  calculating  f  for  any  stage  n ,  where  n*N  ,  is 
determining  the  state  J  at  the  next  stage  (  n+1  ).  Recall,  the  state  at  the 
next  stage  is  calculated  using  the  transition  equation.  Unfortunately,  there 
is  no  standard  form  for  the  transition  equation.  It  is  usually  dependent 
upon  the  specific  problem. 


To  apply  dynamic  programming  to  a  queuing  system,  the  following 
must  be  defined:  the  system  stage  (n),  system  state  (i),  the  stage  decision 
(kn),  the  return  function  (R),  the  recursion  equation  (f),  and  the  transition 
equation.  Generally,  the  system  stage  (n)  is  some  specified  time  period,  e.g. 
n=l  corresponds  to  hour  one,  n=2  corresponds  to  hour  two,  and  so  on. 
The  system  state  in  a  queuing  system  is  usually  the  number  of  customers  in 
line,  and  the  stage  decision  is  the  number  of  servers  to  have  open  for  the 
stage.  Definition  of  the  return  function  is  not  so  simple.  One  way  is  to  use 
a  cost  function  (Magazine,  1971:178): 


K(i)+(ka-kn-i)A+knC 

.KdMkn-i-kjB+knC 


if  kQ>  kn-i 


ifka<kn-i 


(3) 


where 

A  -  cost  of  opening  a  server 
B  =  cost  of  closing  a  server 

C  *  cost  of  operating  an  open  server  for  one  time  period  (stage) 

K(i)  =  holding  cost,  i.e.  cost  incurred  when  i  customers  are 
observed  in  the  system 


A,  B,  and  C  are  usually  fairly  easy  to  determine.  However,  it  is  extremely 
difficult  to  define  the  customer  holding  cost,  K(i).  This  customer  holding 
cost  is  basically  an  attempt  to  put  a  dollar  value  on  the  number  of 
customers  in  line.  The  rationalization  for  this  is  that  if  the  line  is  too  long, 
business  will  be  lost.  The  holding  cost  is  a  measure  of  the  amount  of  the 
lost  business.  Obviously,  any  value  used  for  K(i)  can  only  be  an  estimate. 

The  computational  difficulty  in  solving  a  dynamic  programming 
problem  is  strongly  related  to  the  number  of  possible  states. 
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Unfortunately,  the  addition  of  a  resource  constraint  (e.g.  a  limit  on  the 
total  hours  available  to  be  scheduled)  can  dramatically  increase  the  number 
of  possible  states.  The  reason  for  this  is  the  way  that  resource  constraints 
are  typically  handled.  Usually,  an  extra  state  variable,  corresponding  to 
the  amount  of  resource  remaining,  is  added  (Denardo,  1982:35).  Recall  the 
prior  example,  where  the  system  state  was  given  as  the  number  of 
customers  in  line  (i).  When  the  resource  constraint  is  added,  the  new 
system  state  is  now  a  combination  of  the  number  of  customers  in  line  (i) 
and  the  amount  of  resource  remaining  (y).  If  the  maximum  number  of 
customers  in  line  is  I  and  the  total  available  increments  of  resource  is  Y , 
then  the  total  possible  number  of  states  is  now  approximately  (I  x  Y)  . 

This  increase  in  the  computational  difficulty  of  a  dynamic  programming 
problem  as  the  number  of  state  variables  is  increased  is  known  as  the 
“curse  erf  dimensionality*  (Bellman,  1957: ix). 

To  avoid  the  curse  of  dimensionality  in  resource  allocation  problems, 
an  alternative  to  dynamic  programming  is  available.  This  algorithm  is 
called  the  “maximal  marginal  return*  procedure  (Larson  and  Casti, 
1982:350).  The  maximal  marginal  return  procedure  starts  with  none  of 
the  resource  allocated.  Each  unit  of  the  resource  is  then  added  so  that  it 
maximizes  the  immediate  marginal  return  (or  for  a  minimization  problem, 
each  unit  of  resource  is  added  to  minimize  the  immediate  marginal  return). 
The  procedure  is  complete  when  all  units  of  the  resource  are  allocated. 

The  algorithm  is  simple  and  is  usually  more  efficient  than  dynamic 
programming.  Unfortunately,  a  condition  for  its  use  is  that  the  return 
function  at  each  stage  is  independent  of  the  return  functions  at  all  other 
stages.  This  condition  is  often  violated  in  a  queueing  problem,  where  the 


state  of  the  system  (number  of  customers  in  line)  depends  upon  the  actions 
taken  at  previous  stages. 

The  Fluid  Approximation  for  Queues  (Kieinrock,  1976:56-02) 

Analysis  of  queueing  systems  is  complex  because  it  involves  several 
random  variables;  time  between  customer  arrivals  is  a  random  variable, 
and  the  time  required  to  serve  a  customer  is  a  random  variable.  This 
research  effort  is  an  attempt  to  control  queue  length  by  varying  the 
number  of  servers  to  the  queue.  To  do  this,  the  effect  of  the  number  of 
servers  upon  the  queue  length  must  be  known. 

Queues  are  generally  classified  by  the  distribution  of  customer 
interarrival  times,  the  distribution  of  service  times,  and  the  number  of 
servers  for  the  queue.  For  some  special  distributions  of  customer 
interarrival  times  and  service  times,  the  exact  relationship  between  number 
of  servers  and  queue  length  can  be  derived.  Probably  the  most  well- 
known  is  the  case  where  the  distribution  of  the  customer  interarrival  times 
is  exponential  and  the  distribution  of  the  service  times  is  exponential  (the 
famous  M/M/s  queue).  For  more  general  cases,  when  the  distributions  are 
unknown  or  not  well-behaved,  approximations  must  be  used  to  obtain  a 
relationship  between  the  number  servers  and  the  queue  length.  The  fluid 
approximation  is  one  such  approximation. 

In  any  queuing  system,  the  number  of  customers  is  a  discontinuous 
function  of  time.  This  is  because  the  number  of  customers  can  only  change 
in  integral  units— half  a  customer  does  not  exist.  The  number  of 
customers  in  the  system  at  time  t  can  be  given  as: 

N(t)  =  A(t)-D(t)  (4) 
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where 

A(t)  =  Number  of  customer  arrivals  in  (0,  t) 

D(t)  =  Number  of  customer  departures  in  (0,  t) 

In  Figure  2,  the  relationship  between  A(t),  D(t),  and  N(t)  is  shown 
graphically. 


(Klemrock,  1976:57) 

Figure  2  Number  of  Customers  in  System  as  a  Function  of 

Arrivals  and  Departures 

The  fluid  approximation  to  a  queue  takes  advantage  of  the  fact  that 
when  a  system  is  in  a  heavy  traffic  condition,  the  number  of  customers  can 
be  represented  as  a  continuous  function  of  time  instead  of  a  discontinuous 
function  of  time.  This  is  accomplished  by  using  the  average  number  of 
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customer  arrivals  and  departures  instead  of  the  exact  values,  giving  the 
fluid  approximation  as: 

Nf(t)  -  AOT  -  W  (5) 

where 

A(t)  =  Mean  number  of  customer  arrivals  in  (0,  t) 

D(t)  =  Mean  number  of  customer  departures  in  (0,  t) 

If  the  customer  arrival  rate  as  a  function  of  time  is  given  as  X(t)  and  the 
service  rate  as  a  function  of  time  is  given  as  fi(t),  then: 


t 


> 

II 

> 

3 

+ 

* 

(6) 

t 

D(t)  =  D(0)  +  fn(y)cty 

(7) 

The  fluid  approximation  is  shown  graphically  in  Figure  3. 

There  are  several  limitations  of  the  fluid  approximation  that  should 
be  noted.  The  primary  assumption  of  the  fluid  approximation  is  that  the 
system  is  in  heavy  traffic.  If  the  queue  empties  out  and  servers  become 
idle,  the  approximation  will  no  longer  be  accurate.  Problems  can  also  arise 
at  the  other  extreme— -a  sudden  large  influx  of  customers.  The  fluid 
approximation  tends  to  underestimate  queue  length  when  the  system 
reaches  saturation  in  a  very  short  time  (the  typical  situation  during  a  rush 
hour,  where  X(t)  *  y,(t)  ). 

Although  the  queue  length  may  be  underestimated  for  some 
situations,  the  fluid  approximation  is  still  valid  when  X(t)  *  sji(t)  .  Most 
attempts  to  schedule  servers  to  queues  use  the  steady-state  results  of 


Figure  3  Fluid  Approximation  to  N(t) 

the  simple  M/M/s  model  of  a  queue  (customer  interarrival  times  and 
service  times  are  distributed  exponentially).  One  problem  with  this  model 
is  that  it  is  only  valid  when: 

p  -  -  *  1  (8) 

which  means  that  the  arrival  rate  cannot  exceed  the  overall  service  rate. 

To  overcome  this  limitation,  the  number  of  servers  s  must  be  chosen  so 
that  Eq  (8)  is  obeyed.  This  is  the  approach  used  by  Segal  (1974)  and  by 
Kwan,  Davis,  and  Greenwood  (1988).  However,  there  is  a  more 
fundamental  problem  involved  with  using  the  steady-state  M/M/s  model. 
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If  the  customer  arrival  rate  and  the  overall  service  rate  are  changing  over 
time,  then  the  system  never  reaches  steady-state  (Kwan,  Davis,  and 
Greenwood,  1988:267). 

The  fluid  approximation  is  useful  because  it  is  essentially  a 

deterministic  approximation  to  queue  behavior.  In  other  words,  if  one 
knows  A(t)  and  D(t)  ,  A(t)  and  D(t)  ,  or  X(t)  and  |t(t),  then  Nf(t)  * 

N  (t)  can  be  determined.  With  most  other  approximations,  all  that  can  be 
deduced  is  the  probability  distribution  of  N(t). 

Shift  Scheduling 

The  staff  scheduling  problem  is  basically  a  problem  of  scheduling  a 
limited  resource— people— to  meet  a  set  erf  requirements  at  a  minimum 
cost.  The  general  formulation  of  this  problem  is  (Baker,  1976:157): 


where 

xj 

cj 

bi 

A 

aU 


Min  cx 
s.t.  Ax  2  b 

x  2  0  ,  integer 

the  number  of  workers  for  shift  j 
the  cost  of  assigning  a  worker  to  shift  j 
the  number  of  worker  required  during  period  i 
a  0-1  matrix  with  elements  ay 

1  if  shift  j  works  during  period  i 
0  if  shift  j  does  aot  work  during  period  i 


(9) 


Problem  (9)  is  a  general  integer  programming  problem  and  can  thus  be 
difficult  to  solve.  However,  a  special  case  of  the  above  problem  is  the  cyclic 
scheduling  problem,  where  each  column  of  A  consists  of  consecutive  ones 


and  zeroes.  Such  a  case  might  occur  when  scheduling  a  5-day  work  week 
where  the  days  off  must  be  consecutive.  For  this  case,  A  would  be  given 
as: 


A 


10  0  1111 
I  1001  I  1 
1110  0  11 
11110  0  1 
1  1  1  I  1  0  0 

o  i  it  1  i  o 

0  0  11111 


(10) 


The  consecutive  ones  in  the  A-matrix  of  Eq  (10),  sometimes  referred  to  as 
circular  ones  (Bartholdi,  1881:503),  indicate  that  a  worker  is  available 
continuously  during  all  consecutive  time  periods.  If  this  assumption  is  met, 
algorithms  are  available  to  solve  the  problem  more  efficiently  than  general 
integer  programming  (Bartholdi,  1981:503). 

Segal  considers  the  case  where  work  periods  are  considered  to  be 
hours  in  the  day  instead  of  days  in  the  week  (Segal,  1974).  For  this  case 
the  A-matrix  will  be  linear  instead  of  cyclic.  That  is,  all  ones  in  a  column 
of  A  will  be  adjacent  (the  top  and  bottom  rows  of  A  are  not  considered 
adjacent).  In  Eq  (ll),  an  example  of  such  a  matrix  is  shown  for  eight  and 
five-hour  work  shifts. 

To  solve  the  problem  given  in  (9),  with  A  as  in  Eq  (1 1),  Segal 
converts  the  problem  into  a  network  as  shown  in  Figure  6  (Segal, 
1974:812-815).  In  Segal's  network  formulation,  the  nodes  correspond  to 
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A 


100010000000 
110011000000 
111011100000 
111111110000 
111101111000 
111100111100 
111100011110 
111100001111 
011100000111 
001  10000001  1 
00010000000) 


(11) 


the  transition  from  one  time  period  to  the  next.  The  forward  arcs  from  i 
to  i+1  correspond  to  the  actual  time  periods.  The  backward  arcs  from  m 
to  l  (1  <  m)  represent  the  possible  shifts.  The  arc  parameters  can  be 
defined  as  follows: 

Forward  Arcs: 

ui,i+l  “  upper  capacity  of  arc  *  the  maximum  number  of  workers 

allowed  at  one  time  plus  an  estimate  of  the  number  of  workers 
on  break 

Li, i+1  =  lower  capacity  of  arc  =  bj 

Ci,i+i  =  0 
Backward  Arcs: 

Umji  =  upper  capacity  of  arc  =  the  maximum  number  of  workers 
available  to  work  this  shift 

Lm,l  =  lower  capacity  of  arc  =  the  minimum  number  of  workers  to 
be  assigned  to  work  this  shift 

Cm,l  *  the  unit  cost  of  this  shift 
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This  problem  can  be  solved  efficiently  using  a  network  flow  algorithm.  If 
the  worker  requirements  (bi)  are  all  integers,  then  the  solution  is 
guaranteed  to  be  integer. 


1900 


1800 


1600 


I*  1  (6 


1300 


l-l  (4 


1 100 


1000 


0900 


Figure  4  Segal's  Network  Conversion 


All  of  the  scheduling  algorithms  discussed  thus  far  assume  that  each 
requirement  bi  gives  the  minimum  number  of  workers  needed  during 
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period  i .  This  is  why  the  constraints  in  (9)  are  given  as  inequalities.  If 
the  requirements  are  actually  ideal  numbers  of  workers  required  during 
each  period,  then  the  constraints  in  (9)  should  be  changed  to  equalities. 
However,  if  this  were  the  only  change  made,  then  the  problem  will  often 
turn  out  to  be  infeasible.  Therefore,  an  integer  goal  programming 
formulation  should  be  used  (Koelling  and  Bailey,  1984:302): 


where 

V 

atj 

xj 

bj 


Xj  ^  0  ,  J  =  1,  ... ,  n 

some  achievement  function  to  be  specified 
as  defined  above 
as  defined  above 
as  defined  above 

number  of  workers  below  requirement  for  time  period  i 
number  of  workers  above  requirement  for  time  period  i 


Again,  as  long  as  the  bj  are  all  integers  and  provided  V  is  linear,  the 
solution  to  (10)  is  guaranteed  to  be  integer.  Baker  gives  an  almost  identical 
formulation  (Baker,  1978:181).  The  only  real  difference  is  that  Baker 
defines  the  objective  function  specifically  as: 

m  n  n 

Minimize  2  cj  * j  +  2  <*i  d*  +  2  Pi  (13) 

j-1  i-1  i-1 
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Network  Programming  (Veinott  and  Wagner,  1962:520) 

It  was  no  coincidence  that  Segal  was  able  to  convert  problem  (9)  into 
a  network  problem.  As  early  as  1962,  Veinott  and  Wagner  showed  that 
any  problem  such  as  (9),  where  each  row  of  the  A  matrix  consists  of 
consecutive  zeroes,  followed  by  consecutive  ones,  followed  by  consecutive 
zeroes,  could  be  converted  into  an  equivalent  network  problem.  Suppose 
that  A  is  m  x  n  .  Each  constraint  except  the  first  is  replaced  by  itself 
minus  the  previous  constraint.  The  first  constraint  is  left  intact.  An 
additional  constraint,  equal  to  the  last  constraint  times  -1  ,  is  added.  This 
method  is  illustrated  by  the  following  example  (adapted  from  Veinott  and 
Wagner,  1962:520): 


■1111000000' 

“Xl“ 

X2 

"bf 

0111111000 

• 

b2 

0011011110 

• 

b3 

0001001011 

• 

Lb*J 

L»  J 

-X10- 

(14) 


After  performing  the  transformation,  one  obtains: 


1  1 

-1  0 

0  -1 
0  0 

0  0 


1  t 
0  0 
0  0 
-l  0 
0  -1 


0  0 
t  1 
-1  0 
0  -1 
0  0 


0  0 

1  0 

0  1 

0  -1 
-1  0 


0  0 
0  0 
1  0 
0  1 
-l  -1 


rxi 

x2 


i-xio-* 


-hi 

b2-bi 

bs-b2 

b4-b3 

-  -b*  - 


(15) 


The  equations  of  (15)  have  the  required  structure  to  be  represented  as  a 
network.  Each  column  of  A  contains  a  single  one,  a  single  minus  one, 
and  all  remaining  entries  are  zero,  and  thus  A  can  be  thought  of  as  the 
node-arc  incidence  matrix  of  a  network. 
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The  Out-of- Kilter  Algorithm  (Fulkerson,  1981:18-27) 

An  efficient  algorithm  for  solving  problems  of  the  form: 


Min  cx 

s.t.  Ax  -  b  (18) 

0  £  X  £  U 


where  A  is  a  node-arc  incidence  matrix  of  a  network,  is  the  out-of-kilter 
algorithm.  The  out-of-kilter  algorithm  is  better  than  the  simplex 
algorithm  for  problems  of  this  form  because  it  eliminates  the  need  to  carry 
the  basis  inverse  and  can  thus  reduce  the  computational  burden  of  solving 
the  problem. 

The  dual  of  problem  (18)  can  be  written: 


Max  bx  -  up 
s.t.  *  A  -  f,  £  c 


ft  2  0 

By  defining: 

F(j)  =  “From"  Node  *  the  originating  node  of  arc  j 
TO)  =  "To"  Node  =  the  destination  node  of  arc  J 

’  *f0)  "  *tQ)  ~  CJ 


the  Kuhn- Tucker  conditions  for  optimality  can  be  reduced  to: 


Ax  -  b 


Cj<0  whenxj=0 


for  arc  j  :  | 


cj=G  whenO£Xj£Uj 
Cj>Owhenxj=Uj 


(17) 


(18) 

(19) 
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For  a  given  ( x,  *  ),  arc  j  is  said  to  be  in-kilter  if  (19)  is  satisfied. 
Otherwise,  arc  j  is  said  to  be  out-of-kilter.  If  an  arc  is  out-of-kilter,  the 
kilter  number  is  the  amount  of  flow  required  to  convert  the  arc  to  the  in¬ 
kilter  condition.  For  ( x,  x  )  to  be  a  solution  to  (16),  the  kilter  number  for 
each  arc  must  be  zero  (all  arcs  must  be  in-kilter). 

The  out-of-kilter  algorithm  consists  of  two  phases.  In  the  primal 
phase,  all  dual  variables  are  held  fixed,  and  the  primal  variables  are 
changed  in  an  attempt  to  reduce  the  sum  of  all  kilter  numbers.  In  the  dual 
phase,  the  process  is  reversed.  The  primal  variables  are  held  fixed,  and 
the  dual  variables  are  changed  in  an  attempt  to  reduce  the  sum  of  all  kilter 
numbers. 


Lagrangian  Relaxation  (Fisher,  1981:1-8) 

Many  difficult  integer  programming  problems  can  be  viewed  as  easy 
problems  complicated  by  a  relatively  small  number  of  side  constraints. 

Such  a  problem  can  be  written  as: 

Z  =  Min  cx 

s.t.  Ax  =  b  (20) 

Dx  =  e 

x  i  0  ,  integral 

where  Ax  =  b  are  the  difficult  constraints.  The  Lagrangian  relaxation  is 
formed  as  follows: 

Zq(u)  =  Min  cx  +  u  (Ax  -  b) 


s.t.  Dx  =  e 

x  i  0 ,  integral 


(21) 


where  u  *  (uj,  U&  ...  ,  Um)  is  a  vector  of  Lagrange  multipliers.  For  (20) 
and  (21),  the  following  inequality  will  always  hold: 

Zq(u)  s  cx*  +  u  (Ax*  -  b)  =■  Z  (22) 

In  general,  it  is  not  possible  to  find  Zd(u)  =  Z  .  However,  if  Zd(u)  =  Z  , 
then  by  (22),  x*  is  optimal. 

Obviously,  for  a  fixed  u  ,  (21)  is  easy  to  solve.  Ideally,  u  should 
be  found  so  that  Zd(u)  *  Z  .  The  best  choice  of  u  is  the  solution  to  the 
problem: 

Zd  *  max  Zd(u)  (23) 

u 

If  Zq(u)  were  differentiable  at  all  points,  u  could  be  found  by  setting  the 
gradient  of  Zd(u)  equal  to  0.  Unfortunately,  Zd(u)  is  not  differentiable  at 
all  points,  so  this  will  not  work  here.  An  adaptation  of  the  gradient 
method,  known  as  the  subgradient  method,  has  become  a  popular  approach 
to  selecting  u  .  In  this  method,  a  sequence  {  u*  }  is  generated  starting  at 
a0  =  0  and  using: 

uk+l  »  +  tfc(Axk  -  b)  (24) 


where  xk  is  an  optimal  solution  to  (21)  at  the  previous  iteration  and  t* 
is  a  positive  step  size.  A  common  equation  for  the  step  size  is: 


tk  - 


m 

2 

i-i 


n 


2aijx/-bi 

Lj-l 


(25) 


r?r  ~ 


where  Xk  is  a  scalar  between  0  and  2.  An  empirical  rule  for  X*  is  to  set 
Xo  =  2  and  set  X*  =  (Xk-i/2)  if  Zd(u)  has  failed  to  increase  in  a  specified 
number  of  iterations.  It  should  be  noted  that  while  Eq  (25)  has  worked 
often  in  practice,  there  is  no  guarantee  that  it  will  always  force 
convergence  to  the  optimal  solution. 

Summary 

This  chapter  presented  a  review  of  literature  for  six  areas  relevant 
to  this  research.  First,  dynamic  programming  with  resource  constraints 
was  discussed.  Second,  the  fluid  approximation  for  obtaining  an  estimate 
of  queue  length  was  reviewed.  The  third  area  of  interest  in  the  literature 
review  were  previous  efforts  at  shift  scheduling.  The  fourth  part  of  the 
review  showed  how  certain  integer  programming  problems  could  be 
converted  into  network  problems,  and  the  fifth  part  briefly  discussed  an 
efficient  algorithm  for  solving  these  network  problems.  The  sixth  and  final 
section  of  the  literature  review  discussed  Lagrangian  relaxation,  a  method 
for  removing  or  relaxing  certain  constraints  that  change  an  otherwise 
easily  solved  problem  into  a  more  difficult  problem. 
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The  scheduling  of  servers  at  a  service  organization  is  a  three  phase 
problem.  In  Phase  I,  the  server  requirements  for  each  time  period  must  be 
determined  so  that  the  desired  customer  waiting  time  is  obtained.  Then  in 
Phase  II,  the  optimal  number  of  workers  to  schedule  to  each  possible  shift 
must  be  found  so  that  the  server  requirements  determined  in  Phase  1  are 
met.  Finally,  in  Phase  III,  actual  workers  are  assigned  to  the  shifts  in  the 
numbers  calculated  in  Phase  II. 

The  Phase  I  problem  is  solved  here  using  dynamic  programming. 

The  output  of  the  Phase  I  dynamic  program  becomes  the  input  to  Phase  II. 
Because  of  the  special  structure  of  the  constraints  on  the  possible  worker 
shifts,  the  Phase  II  problem  can  be  solved  using  a  network  flow  algorithm. 
The  Phase  III  problem  was  not  addressed  here.  One  approach  to  the  Phase 
III  problem  might  be  to  view  it  as  an  assignment  problem,  where  workers 
are  assigned  to  the  shifts  calculated  in  Phase  II. 

The  solution  to  the  Phase  I  and  Phase  II  problems  were  implemented 
using  Turbo  Pascal.  The  resulting  code,  named  the  Server  Scheduler,  is 
given  in  Appendix  A. 

Phase  I 

The  purpose  of  Phase  I  is  to  determine  the  checker  requirements 
throughout  the  month  to  obtain  a  mean  customer  waiting  time  of  five 
minutes.  The  total  checker  hours  must  be  less  than  a  pre-set  number  of 
hours.  Customer  arrivals  to  the  queue  have  been  tabulated  throughout  Air 
Force  commissaries  at  one-hour  intervals.  Because  of  the  time-staged 
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nature  of  the  arrival  data,  a  dynamic  programming  approach  is  suggested 
for  determining  the  checker  requirements  throughout  the  month.  The 
dynamic  programming  approach  assumes  that  the  customer  arrival  rate 
and  the  service  rate  can  be  treated  deterministically  (the  validity  of  this 
assumption  is  discussed  in  Chapter  IV).  It  has  been  shown  that  the  service 
time  for  a  customer  is  given  by  (Moulder,  1987:105): 

service  time  *  1.48  +  0.05  y  (28) 

where  y  is  a  random  variable  having  a  gamma  distribution  with 
parameters  a  =  3.2  and  0  =  23.1.  The  mean  value  of  y  is  73.92,  so  the 
mean  service  time  is  5. 158  minutes.  The  corresponding  service  rate  fi  is 
0. 194  customers  per  minute. 

The  natural  stages  in  this  problem  are  each  hourly  interval  (n).  It 
was  not  possible  to  formulate  a  return  function  in  such  a  way  that  the 
return  functions  at  each  stage  would  be  independent,  and  so  the  faster 
maximal  marginal  return  method  could  not  be  used.  If,  however, 
conventional  dynamic  programming  causes  the  maximum  total  checker 
hours  to  be  exceeded,  then  a  variant  of  the  marginal  return  method  will  be 
used  to  deallocate  checkers.  As  in  most  queueing  problems,  the  state 
variable  is  chosen  as  the  number  of  customers  in  line  (Ln).  The  decision 
variable  is  the  number  of  checkers  (kn)  to  have  open  during  period  n.  As 
mentioned  in  Chapter  11,  it  is  usually  difficult  to  formulate  a  good  return 
function  for  a  queueing  problem.  This  is  not  the  case  here.  Because  the 
Air  Force  Commissary  Service  has  requested  that  the  mean  customer 
waiting  time  be  kept  at  five  minutes  throughout  the  day,  the  return 
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function  can  be  defined  as  the  deviation  of  the  waiting  time  from  five 
minutes,  i.e. 


Rn  =  lWn-5  I 


(27) 


where  Wn  is  the  mean  customer  waiting  time  during  period  n.  The 
mean  customer  waiting  time  can  be  defined  as: 


Wn 


i-n+  Ln+l 
2likn 


(28) 


where 

Ln  =  number  of  customers  in  line  at  stage  n 
Ln+i  =  number  of  customers  in  line  at  stage  n+1 
)i  =  mean  service  rate 
kn  *  number  of  checkers  at  stage  n 


In  Eq  (27),  [(Ln  +  Ln+i)  /  2]  gives  the  average  number  of  customers  in  the 
one-hour  interval  while  }ikn  gives  the  overall  service  rate  of  all  checkers 
combined.  The  number  of  customers  in  line  at  stage  n+1  is  related  to  the 
number  of  customers  in  line  at  stage  n  by. 


Ln+i 


L„  +  Aft  -  60)ikn  if  n  *  final  hour  of  day 
0  if  n  *  final  hour  of  day 


(29) 


where  all  variables  are  as  defined  above,  and  An  is  the  number  of 
customer  arrivals  during  the  one-hour  Interval  (which  has  been  tabulated). 
Eq  (29)  is  a  fluid  approximation  to  the  behavior  of  the  queue,  with  arrivals 
A(t)  =  An  and  departures  D(t)  *  60)ikn  .  The  case  of  Ln+i  =  0  is  used 
to  force  the  queue  to  empty  at  the  end  of  each  day  and  start  empty  for  the 
following  day.  For  this  to  occur,  the  departures  in  the  last  period  of  the 


day  must  exceed  the  total  number  of  leftover  customers  plus  the  number  of 

i 

arrivals  in  the  last  period,  i.e.  60}ikn+i  *  Ln  +  An+i  .  Thus  the  normal 
flow  of  a  dynamic  program  as  shown  in  Figure  1  is  modified  as  shown  in 
Figure  5. 

The  only  other  formula  needed  to  complete  the  dynamic 
programming  formulation  is  the  backward  recursion  formula.  It  is  as 
given  in  Eq  (l),  where  the  objective  is  minimization: 


fn  * 


min  £  Rn+  fn+1  ]  for  n<  N 

min  [RnI 
*N  L  1 


for  n=N 


(30) 


The  dynamic  programming  formulation  of  the  Phase  I  problem  can 
be  summarized  as  follows: 


Stage  Variable 

n  =  each  one-hour  interval 
State  Variabte 

Ln  =  number  of  customers  in  line 
Qseision  Variable 

kn  =  number  of  checkers  open 

Transition  Equation 


An  -  60|tkn  if  a  *  final  hour  of  day 
0  if  n  «  final  hour  of  day 

and  60lLkn  a  Ln*|  ♦  An 


Return  Function 

Rn  *  lWn-5| 


Ln+Ln+1 

2likn 

2Ln~t~An-60jtkn 


2ltkn 


-  5 
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[ftn+fn+l  ] 

M 


for  n<N 
for  n=N 


The  Phase  1  dynamic  programming  technique  is  implemented  in  the 
Server  Scheduler  in  the  procedure  “Phasel.*  This  procedure  in  turn  makes 
calls  to  procedures  to  calculate  Ln+i  (<rcalc_num_custJ,)}  Wn 
(“waiting_time,0,  and  Rn  (“Return*). 

The  output  of  Phase  I  is  a  the  vector  k  containing  the  checker 
requirements  needed  throughout  the  month  to  achieve  a  five  minute 
waiting  time.  If  the  total  checker  hours  exceed  the  maximum  total  hours, 
i.e. 

N 

2  kn  >  Max  Total  Hours  (31) 

n-l 


then  checkers  must  be  removed  until  the  total  checker  hours  are  equal  to 
the  maximum  total  hours.  A  variation  of  the  maximal  marginal  return 
method  is  used  to  determine  the  stage  from  which  to  remove  a  checker. 
Checkers  are  removed,  one  at  a  time,  from  the  stage  that  produces  the  least 
gain  in  the  overall  sum  of  the  return  functions.  In  the  Server  Scheduler, 
checker  removal  via  the  marginal  return  method  is  accomplished  in  the 
procedure  “Deallocate.  * 
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Phase  11 

Once  the  checker  requirements  are  calculated  in  Phase  I,  the  optimal 
shift  schedules  can  be  determined.  Because  of  the  number  of  possible  shifts, 
it  is  best  to  solve  Phase  II  for  each  day  separately  and  then  combine  the 
daily  shift  schedules  into  an  overall  monthly  schedule.  It  is  assumed  that 
possible  shift  lengths  are  4,  5,  6,  7,  or  8  hours  long. 

The  Phase  II  problem  is  formulated  as  in  equations  (12)  and  (13): 


Min  ld+  +  ld" 

s.t.  Ax  +  d~  -  d+  »  b  (32) 

x  i  0  ,  integer 
d“  i  0  ,  integer 
d+  i  0  ,  integer 


where  all  variables  are  as  defined  in  Chapter  II.  A  sample  constraint 


matrix  for  problem  (32)  is  given  in  Eq  (33).  For  simplicity,  the  sample 


constraint  matrix  shows  only  shift  lengths  of  4  and  8  hours. 
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(33) 


Using  the  method  of  Veinott  and  Wagner  (1962)  as  outlined  in  Chapter  II, 


the  constraints  are  transformed  into  a  network  node-arc  incidence  matrix. 
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The  constraints  are  illustrated  in  Eq  (34),  and  the  corresponding  network 
is  shown  in  Figure  6. 
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Since  the  network  problem  given  in  Eq  (34)  and  Figure  6  is  equivalent  to 
the  problem  given  in  Eq  (33),  the  solution  to  the  network  problem  is  the 
solution  to  the  original  integer  goal  programming  problem.  But  the 
solution  to  the  network  problem  can  be  found  much  more  quickly  than  the 
solution  to  the  integer  programming  problem.  More  importantly,  since  the 
elements  of  b  are  all  integer,  the  solution  to  (32)  is  guaranteed  to  be 
integral  if  the  constraints  are  as  in  Eqs  (33)  and  (34).  In  the  Server 
Scheduler,  the  out-of-kilter  algorithm  is  used  to  solve  the  equivalent 
network  programs.  The  procedure  used  is  called  “MinCostFlow,*  and  is 
based  on  the  implementation  of  the  out-of- kilter  algorithm  given  in 
Kennington  and  Helgason  (1980:78-88). 
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Lunch  Breaks 

Because  a  special  network  flow  algorithm  can  be  used  to  solve  it,  the 
Phase  II  procedure  outlined  above  is  very  efficient.  However,  this 
formulation  has  not  made  provision  for  workers  to  take  lunch  breaks. 
Assuming  that  a  worker  on  any  shift  longer  than  six  hours  is  entitled  to  a 
one-hour  lunch  break,  the  sample  constraint  matrix  shown  in  Eq  (33)  must 
be  modified  as  shown  in  Eq  (35): 
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(35) 


This  new  constraint  matrix  is  no  longer  amenable  to  the  network 
transformation,  and  the  solution  to  the  problem  is  no  longer  guaranteed  to 
be  integral.  Therefore,  to  solve  this  problem  one  must  resort  to  more 
difficult  general  integer  programming  techniques.  For  a  typical  day,  when 
five-,  six-,  and  seven-hour  shifts  are  added,  the  problem  given  in  (35) 
would  have  50  to  60  integer  variables  that  range  from  0  to  30.  Each 
individual  variable  would  need  five  0-1  variables  (25=32).  This  gives  a 
total  of  250  to  300  variables  in  a  0-1  integer  programming  formulation 
(too  large  to  solve  on  a  microcomputer). 

Several  alternate  formulations  are  available  that  can  help  avoid  the 
problems  associated  the  formulation  of  (35).  The  first  alternate 
formulation  adds  additional  constraints  for  every  hour  corresponding  to  a 
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shift  lunch  hour.  These  constraints  are  then  dualized  to  form  a 
Lagrangian  relaxation  to  the  original  problem.  The  Lagrangian  relaxation 
retains  the  consecutive  ones  in  the  constraint  matrix,  thus  guaranteeing  an 
all-integer  solution.  The  second  alternate  formulation  splits  each  eight- 
hour  shift  into  two  four- hour  shifts,  each  seven-hour  shift  into  a  four  and 
a  three  hour  shift,  and  each  six-hour  shift  into  two  three  hour  shifts. 

This  formulation  also  retains  the  consecutive  ones  in  the  constraint  matrix, 
and  the  solution  is  again  guaranteed  to  be  all- integer. 

Lagrangian  Relaxation  Formulation.  In  this  formulation,  as  in  the 
formulation  of  (35),  the  length  of  a  shift  with  a  lunch  hour  is  increased  by 
one  hour.  An  eight-hour  shift  is  changed  into  a  nine-hour  shift,  a  seven- 
hour  shift  is  increased  into  an  eight-hour  shift,  and  a  six-hour  shift 
becomes  a  seven-hour  shift.  Unlike  the  formulation  of  (35),  the  added  hour 
in  a  shift  is  not  a  zero  in  the  constraint  matrix.  The  new  formulation  is 
given  in  (36): 
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The  constraint  matrix  of  (38)  is  not  complete.  If  the  only  the 
constraints  of  (36)  were  used,  the  number  of  checkers  during  the  fifth, 
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sixth,  and  seventh  hours  would  be  overestimated  by  xi,  and  X3 
respectively.  The  reason  for  this  is  that  the  checkers  assigned  to  these 
shifts  are  actually  out-to-lunch  during  these  hours.  This  is  best  illustrated 
by  example.  During  the  fifth  hour,  the  constraints  of  (36)  indicate  that  the 
number  of  checkers  on  duty  is: 

X1  +  X2  +  X3  +  X5  +  X8  +  X7  +  X8  (37) 

Since  checkers  from  the  first  shift  are  out  to  lunch,  the  actual  number  of 
checkers  on  duty  is: 

X2  +  X3  +  X5  +  X8  +  X7  +  X8  (38) 

In  other  words,  requirement  bs  will  be  undershot  by  xi  .  To  off-set  the 
overestimation  of  checkers,  extra  constraints  are  added  to  force  extra 
checkers  to  be  scheduled  during  lunch  hours.  For  the  example  given 
above,  the  extra  constraint  is  of  the  form: 

■  xi  or  dj  -  xi  =  0  (39) 

The  effect  of  the  added  constraint  is  to  force  xi  extra  checkers  to  be 
scheduled  during  the  fifth  hour.  These  extra  checkers  exactly  offset  the 
shortage  created  by  the  checkers  that  are  out  to  lunch.  The  complete 
formulation  for  the  sample  problem  is  given  in  (40): 
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(40) 
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The  extra  lunch-hour  constraints  in  (40)  destroy  the  consecutive  ones 
property  that  allowed  the  problem  to  be  converted  into  a  network  (and 
insured  an  integer  solution).  To  regain  the  consecutive  ones  property,  the 
lunch-hour  constraints  are  dualized  and  a  Lagrangian  relaxation  is 
formed: 
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Theoretically,  the  problem  formulated  in  (41)  should  be  solvable  by 


the  method  outlined  in  Chapter  2.  First,  the  multipliers  are  set  at  u°  - 
(0,0,0) ,  and  the  resulting  problem  can  be  solved  using  the  network  flow 


algorithms  already  outlined.  Then  the  multipliers  are  updated  using  Eqs 
(24)  and  (25)  and  the  process  is  repeated  until  a  solution  is  found  that 
satisfies  the  extra  lunch  break  constraints  ( tk  -  0  in  Eq  (25) ). 


Unfortunately,  implementation  of  the  Lagrangian  relaxation  method  for 
this  problem  proved  to  be  difficult.  Apparently,  this  is  one  application  for 
which  Eq  (25)  failed  to  produce  convergence  to  the  optimal  solution. 
Lagrangian  relaxation  is  still  a  promising  method  for  solving  the  shift 
scheduling  problem  with  additional  constraints  (including  but  not  limited  to 
lunch  break  constraints).  However,  a  different  policy  must  be  found  to 
update  the  multipliers,  since  Eq  (25)  has  proven  to  be  ineffective  for  this 


application. 


Split-Shift  Formulation.  Another  promising  way  to  account  for 
lunch  breaks  is  to  break  each  shift  requiring  a  lunch  break  into  two  shifts. 
For  example,  an  eight-hour  shift  extending  from  0900  to  1800  with  a  lunch 
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break  at  1300  would  become  two  four-hour  shifts,  the  first  extending  from 
0900  to  1300  and  the  second  from  1400  to  1800.  Similarly,  seven-hour 
shifts  are  split  into  a  four  and  a  three-hour  shift,  and  six-hour  shifts 
become  two  three-hour  shifts.  Using  this  formulation,  the  final  constraint 
matrix  consists  of  five,  four,  and  three-hour  shifts: 
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\  * 

00001  1 1 000001 1  1 0000001 1  I 00000 00000000000  l-l  0000 

b9 

000001  1 0000001  1 00000001  1 0000000000000000001-1  00 

bio 

\0000001  OOOOOOOI 000000001 00000 000 000000 000 000 l-l  y 

\bl  1/ 

The  split-shift  formulation  is  not  without  complications.  First, 
solution  of  the  problem  given  in  (42)  gives  optimal  numbers  of  five-,  four-, 
and  three-hours  shifts.  Some  way  must  be  found  to  convert  these  shifts 
back  into  eight-,  seven-,  six-,  five-,  and  four-hour  shifts.  Since  this 
conversion  could  be  done  by  hand  if  need  be,  this  is  not  a  serious  limitation. 
But  a  serious  limitation  does  exist.  If  three-hour  shifts  are  not  allowable, 
one  must  ensure  that  all  three-hour  shifts  in  the  optimal  solution  can  be 
converted  into  six,  seven,  and  eight-hour  shifts.  There  are  at  least  two 
possible  ways  to  ensure  this.  The  first  way  maintains  the  network 
structure  but  does  not  guarantee  that  a  solution  could  be  found,  while  the 
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second  guarantees  a  solution  but  requires  a  Lagrangian  relaxation  problem 
to  be  solved. 


The  first  method  for  ensuring  the  conversion  of  three-hour  shifts  is 
similar  to  branch-and-bound  methods  for  solving  integer  programs.  The 
problem  of  (42)  is  first  solved  as  given  using  the  network  algorithm  (as 
before).  If  all  three-hour  shifts  can  be  converted  into  six-,  seven-,  and 
eight-hour  shifts,  the  problem  is  solved.  If,  however,  there  are  M3  extra 
three-hour  shifts  that  can  not  be  converted  from  shift  Xm  ,  then  an  upper 
bound  of  X3  -  M3  is  placed  on  X3 ,  and  the  problem  is  solved  again.  The 
process  is  continued  until  a  feasible  solution  is  obtained. 

The  second  method  for  ensuring  that  the  three-hour  shifts  can  be 
converted  consists  of  adding  extra  constraints.  For  each  three-hour  shift, 
all  other  shifts  that  could  combine  with  the  three-hour  shift  to  form  a 
six-,  seven-,  or  eight-hour  shift  are  enumerated.  The  sum  of  the  checkers 
assigned  to  these  other  shifts  must  exceed  the  number  of  checkers  assigned 
to  work  the  three-hour  shifts.  For  the  problem  of  (42),  this  means  that 
the  following  constraints  must  be  added: 


X16  S  X12  +  X20 

(43) 

X17  *  X13  +  X21 

(44) 

X18  *  X14  +  X22 

(45) 

X19  *  X15  +  X23 

(46) 

X20  *  XJ6  +  X24 

(47) 

X21  *  X8  +  X17 

(48) 

X22  *  X9  +  X18 

(49) 

X23  *  Xio  +  X19 

(50) 
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X24  *  XU  +  X20 


(51) 


These  extra  constraints  again  destroy  the  network  structure  of  the 
problem.  To  regain  that  structure,  one  can  again  apply  a  Lagrangian 
relaxation,  with  the  same  complications  experienced  in  the  original 
Lagrangian  relaxation  formulation.  Of  course,  if  three  hour  shifts  are 
allowable,  then  the  formulation  given  in  (42)  can  be  used,  and  the  optimal 
solution  would  be  obtained. 

Summary 

The  Server  Scheduler  Program  implements  a  two-phase  algorithm 
for  scheduling  servers  at  a  service  organization  (specifically,  checkers  at 
U.S.  Air  Force  commissaries).  In  Phase  I  of  the  algorithm,  dynamic 
programming  is  used  to  find  the  number  of  servers  required  during  each 
scheduling  period  to  obtain  a  target  customer  waiting  time.  A  fluid 
approximation  to  queue  length  is  used  to  calculate  the  average  customer 
waiting  time  during  each  period.  Then,  in  Phase  II  of  the  algorithm,  the 
optimal  number  of  servers  to  schedule  to  each  possible  shift  is  found  so 
that  the  requirements  determined  in  Phase  I  are  met.  Integer 
programming  was  used  to  implement  Phase  II  of  the  scheduling  algorithm. 
Because  of  the  special  structure  of  the  constraints  in  the  Phase  II  integer 
program,  network  techniques  could  be  used  to  solve  Phase  II  efficiently. 

Scheduling  of  checker  lunch  breaks  is  a  difficult  extension  to  the 
Phase  II  shift  scheduling  problem.  A  brute  force  approach,  simply  adding 
zeroes  to  the  corresponding  row  of  each  shift,  destroys  the  structure  of  the 
problem.  With  this  approach,  applying  linear  programming  to  the  problem 
is  not  guaranteed  to  produce  an  integer  solution,  so  more  difficult  and 
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time-consuming  integer  programming  methods  must  be  used.  Use  of 
Lagrangian  relaxation  regains  the  special  structure  that  allows  efficient 
network  techniques  to  be  used  iteratively  to  solve  the  problem.  However, 
commonly  used  multiplier  update  formulas  do  not  work  for  this  problem, 
so  more  research  is  needed  to  determine  if  a  Lagrangian  relaxation 
technique  can  be  successfully  applied  here.  Finally,  a  split-shift 
formulation  was  explored.  In  this  formulation,  each  shift  with  a  lunch 
break  was  split  into  two  shifts.  Again,  the  special  network  structure  was 
regained  that  would  allow  this  formulation  to  be  solved  efficiently.  If 
three-hour  shifts  are  allowable,  this  formulation  works.  However,  if 
three-hour  shifts  are  not  allowed,  then  Lagrangian  relaxation  techniques 
must  be  used  to  solve  this  formulation. 
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This  chapter  is  concerned  with  the  validity  of  the  Phase  I  part  of  the 
Server  Scheduler  model,  i.e.  the  determination  of  the  optimal  checker 
requirements.  Two  different  types  of  model  validation  are  discussed  here. 
First,  face  validity  of  the  model  is  briefly  explored.  That  is,  does  the  model 
and  its  output  seem  to  make  sense?  After  face  validity  is  checked,  a 
simulation  is  used  to  see  if  the  model  achieves  its  goals— specifically  the 
achievement  of  the  desired  mean  customer  waiting  time.  The  simulation 
can  also  be  used  to  investigate  the  general  behavior  of  the  system  under 
some  typical  conditions. 

The  data  used  to  test  the  model  was  tabulated  for  three  weeks  in 
May  and  June  of  1987  at  the  Lackland  AFB  Commissary.  It  consisted  of 
hourly  counts  of  the  number  of  customers  arriving  to  the  queue  and  the 
number  of  customers  currently  in  the  queue.  This  data  is  given  in 
Appendix  B. 


Using  the  three  weeks  of  data  given  in  Appendix  B  for  customer 
arrivals,  Phase  I  of  the  Server  Scheduler  was  run  to  determine  the  optimal 
checker  requirements.  In  Figure  16,  the  checker  requirements  kn  and  the 
arrivals  An  are  plotted  against  n  for  a  typical  day.  As  would  be 
expected,  when  the  number  of  customer  arrivals  increases,  more  checkers 
are  required.  Also,  when  the  number  of  customer  arrivals  in  a  period  is 
very  large,  extra  checkers  are  scheduled  in  the  preceding  periods  in  an 
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attempt  to  prepare  for  the  surge.  So  from  a  simple  face  validity 
standpoint,  the  output  of  the  Server  Scheduler  is  consistent  and  sensible. 


Figure  7  Checker  Requirements  for  11  Jun  87 
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To  truly  test  the  performance  of  the  Server  Scheduler,  simulations 
were  constructed  for  two  typical  days  from  the  three- week  period.  The 
days  chosen  were  6  Jun  87  and  11  Jun  87.  Four  different  simulations 
were  run  for  each  day,  each  using  different  assumptions.  A  summary  of 
the  assumptions  for  each  of  the  five  simulations  follows: 

1.  All  parameters  were  assumed  to  be  deterministic. 

Customers  were  assumed  to  arrive  at  a  constant  rate 
throughout  each  hour,  and  the  service  time  for  each 
customer  was  assumed  to  be  5. 156  minutes. 

2.  Customer  interarrival  times  are  distributed 
exponentially  with  a  known  mean.  Service  times  are 
also  randomly  distributed,  with  the  distribution  given  by 
Eq  (24). 

3.  Customer  interarrival  times  are  distributed 
exponentially,  but  now  the  mean  interarrival  time  is  a 
random  variable  with  a  normal  distribution.  The  mean 
of  the  mean  interarrival  time  is  known,  and  the 
standard  deviation  of  the  mean  interarrival  time  is  5%  of 
the  mean.  Service  time  distributions  are  still  given  by 
Eq  (24). 

4.  Same  as  case  3,  but  now  the  standard  deviation  of  the 
mean  interarrival  time  is  10%  of  the  mean. 

The  simulations  were  run  using  SLAM  on  a  DEC  VAX-8650  running 
under  the  VMS  operating  system.  The  SLAM  model  is  shown  in  Figure  8, 
and  the  SLAM  code  and  output  for  each  of  the  8  runs  is  given  in  Appendix 
C.  The  customer  waiting  time  during  each  period  was  averaged  across  50 
runs  of  each  simulation.  A  plot  of  this  average  is  given  in 


Figure  8  SLAM  Model 


Figures  11  to  18  for  each  simulation,  along  with  the  expected  results 
according  to  the  fluid  approximation  used  in  the  Server  Scheduler  (Figures 
9  and  10). 

Comparison  of  the  simulation  results  against  the  expected  results 
according  to  the  fluid  approximation  of  the  Server  Scheduler  shows  that 
for  some  cases,  the  Server  Scheduler  is  fairly  accurate,  and  for  other  cases, 
the  Server  Scheduler  is  very  inaccurate.  Specifically,  the  Server  Scheduler 
produces  good  results  early  in  the  day.  However,  toward  the  end  of  each 
day,  the  customer  waiting  time  in  the  simulations  exhibited  a  marked 
increase.  Two  approximations  were  made  in  the  Server  Scheduler  that 
might  have  caused  this  problem.  They  were  the  approximations  used  to 
calculate  the  number  of  customers  at  the  end  of  each  period  and  to  calculate 
the  mean  customer  waiting  time  for  each  period. 


Average  Cust  omerWaitin  q  Tme 


Figure  9  Fluid  Approximation  for  Entire  Scheduling  Period 
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Figure  10  Fluid  Approximation  for  6  and  11  Jun  87 
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Stochastic  Mean  Arrival  Rati 
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Figure  13  Case  3:  6  Jun  87 
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Figure  14  Case  4:  6  Jun  87 
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!1  Jun  8? 
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Figure  16  Case  2:  11  Jun  87 
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Figure  17  Case  3:  11  Jun  87 


H  Jun  67 

Stochastic  Mean  Arrival  Rate 
(Standard  Deviation  *  103) 


Figure  18  Case  4:  11  Jun  87 
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The  approximation  for  the  mean  customer  waiting  time  of  each 
period  is  given  by  Eq  (28): 


Wn 


Ln+  Lq+1 


(28) 


To  test  the  accuracy  of  this  approximation,  the  Case  2  simulation  for  11 
Jun  87  was  repeated,  and  the  number  of  customers  at  the  end  of  each 
period  was  averaged  across  50  simulation  runs.  This  allows  a  comparison 
between  the  mean  customer  waiting  times  of  each  period  and  the 
approximation  of  Eq  (28)  using  actual  simulation  values  for  Ln  and  Ln+i . 
This  comparison  is  made  in  Figure  19,  and  the  approximation  appears  to 
be  fairly  accurate. 


Figure  19  Comparison  of  Approximation  and 
Actual  Customer  Waiting  Times 


The  approximation  for  number  of  customers  in  line  at  the  end  of 
each  period  is  given  by  Eq  (29): 

_  Ln+An-60|ikn  if  n*finalhour  of day  . 

La+1  ”  0  ifn=finalhourofday  ^ 

The  values  of  Ln  from  the  last  simulation  are  compared  to  the  expected 
values  of  Ln  given  by  the  approximation  of  Eq  (29)  in  Figure  20.  Notice 
that  near  the  end  of  the  day  the  approximation  for  Ln  becomes 
inaccurate.  This  is  the  probable  cause  of  the  long  waiting  times  observed 
at  the  end  of  each  day  in  the  simulations. 


Figure  20  Comparison  of  Approximation  and 
Actual  Customer  Line  Lengths 


To  correct  for  the  end  of  the  day  underestimation  of  line  lengths 
show  in  Figure  20,  an  additional  checker  can  be  added  at  period  10.  The 
effect  of  the  additional  checker,  shown  in  Figure  21,  is  dramatic.  The 
mean  customer  waiting  time  is  reduced  to  acceptable  levels. 


Figure  2t  Effect  of  Additional  Server  at  End  of  Day 


54 


The  primary  objective  of  this  research  was  to  develop  an  analytical 
model  that  could  be  used  to  optimally  schedule  commissary  checkers  so  that 
the  expected  customer  waiting  time  remains  relatively  constant  throughout 
the  day.  Such  a  model  was  developed  and  proved  valid  for  most  of  the 
day.  There  were  slight  problems  with  the  model  at  the  end  of  each  day, 
but  these  could  be  corrected  by  heuristically  adding  a  server  when  the 
customer  waiting  time  begins  to  increase.  Since  commissary  managers 
typically  have  a  number  of  discretionary  employees  available  for  temporary 
surges,  this  model  could  work  in  practice. 

A  subobjective  was  to  make  the  model  sufficiently  general  to  be  used 
in  any  Air  Force  commissary.  The  model  developed  is  actually  general 
enough  to  be  used  in  any  service  organization  where  the  customer  arrival 
pattern  is  known.  It  allows  the  user  to  specify  a  mean  customer  waiting 
time  goal,  the  second  subobjective  of  the  research.  The  deallocation 
procedure  ensures  that  the  third  subobjective,  to  keep  the  total  checker- 
hours  below  a  maximum  allowable  number,  is  met. 

The  checker  scheduling  model  was  validated  using  simulation.  The 
simulation  showed  that  the  model  worked  for  most  of  the  scheduling 
period.  However,  as  mentioned  previously,  the  model  did  tend  to  produce 
overly  optimistic  estimates  of  mean  customer  waiting  times  late  in  the  day. 

Recommendations 


The  Server  Scheduler  program  could  be  a  valuable  tool  for 
scheduling  commissary  checkers.  The  program  is  overly  optimistic  toward 


the  end  of  each  day,  but  this  problem  can  be  corrected  if  only  a  single 
discretionary  employee  is  available  for  one  scheduling  period  near  the  end 
of  the  day.  To  eliminate  the  need  for  discretionary  employees,  a  different 
approximation  could  be  investigated  that  can  successfully  predict  customer 
line  lengths  as  they  approach  zero. 

The  Server  Scheduler  program  implements  two  phases  of  a  three 
phase  problem.  Although  all  objectives  of  the  research  were  met,  there  are 
two  obvious  extensions  of  this  work.  The  first  is  to  solve  the  third  phase 
of  the  problem,  which  is  the  assignment  of  actual  people  to  specific  shifts. 
The  second  is  to  improve  the  Phase  II  part  of  the  algorithm  (where  the 
number  of  checkers  to  assign  to  each  shift  is  determined).  Currently,  the 
Server  Scheduler  does  not  take  lunch  breaks  into  account.  Several  possible 
techniques  for  accomplishing  this  were  outlined  in  Chapter  III,  along  with 
limitations  of  these  techniques. 

Conclusions 

The  Server  Scheduler  provides  a  reliable  and  automated  method  for 
scheduling  commissary  checkers.  The  method  can  save  management  time, 
and  it  can  save  the  commissary  money  through  improved  utilization  of 
checkers.  The  method  is  sufficiently  versatile  to  be  used  in  any  commissary 
in  the  Air  Force. 
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Appendix  A:  Server  Scheduler 


progroa  Server-Scheduler; 


Phase  I  of  this  prograa  calculates  the  optical  nuaber  of 
checkers  needed  throughout  the  day  (eonth)  in  order  to  achieve 
a  five  einute  customer  eaiting  tiee. 


Phase  1 1  of  this  prograe  then  detereines  the  optieal  number  of 
checkers  to  assign  to  each  shift  in  order  to  eeet  the  Phase  I 
checker  requirements. 


uses  PasPr inter; 


const 

MHXSTHQES  ■  372; 
hU  «  0.10395; 

MU  I  MU  -  5. 156; 

TN  «  60.0, 

MRXjCHECKERS  «  30; 
MINjCHECKERS  «  1; 

LMflX  ■  100; 

IMF IMITV  -  3.4E38; 
IMT-INFIMITY  -  2000000000; 
Maxflrcs  «  255;  {  eaxi 

MaxNodes  ■  254;  {  eaxi 


{  eax  •  stages  in  a  eonth  > 
{  average  service  rate  > 
{  average  service  tiee  > 
{  tiee  in  one  stage  } 
nueber  of  open  checkers  } 
nueber  of  open  checkers  } 
{  eax I bub  nuaber  of  customers  in  line  > 
<  biggest  atloeed  real  nuaber  *  infinity  } 
(  biggest  allowed  integer  ■  infinity  } 
*  arcs  allowed  in  HinCostFlow  procedure  } 
*  nodes  allowed  in  HinCostFlow  procedure  } 


I 

lint 


type 

Intflrrayl  »  array 10. .LMflX]  of  Integer; 

Rea I Array 1  ■  arraytO. .LMflX]  of  real; 

Intflrray2  -  array M . .MflXSTflGESl  of  integer; 

DayOfUeek  ■  (Monday, Tuesday, Wednesday, Thursday, Friday, Saturday, Sunday ); 
Hours  *  record 


open  :  integer; 
close  :  integer; 
end; 

date  *  record 
month  :  1 . . 12; 
day  :  1. .31; 
year  :  I960.. 2100; 
end; 

flrcflrray  ■  array 11. .Maxflrcs]  of  integer; 

Modeflrroy  >  array M. .MaxNodes)  of  integer; 


var 

n  :  integer; 
hour  :  integer; 

MM  :  integer; 
TotalStages  :  integer; 
TotalHours  :  integer; 
MaxHours  :  longint; 
k  :  Intflrray2; 
d  :  array! 1 . .MflXSTflGESl 
f  :  array Cl. .MflXSTflGESl 


{  stage  variable  } 
{  each  hourly  interval  in  a  day  } 
{  total  •  stages  in  a  eonth  } 
{  also  total  •  stages  in  a  eonth  > 
{  total  •  checker  hours  allocated  } 
{  eaxi bub  »  checker  hours  to  be  al located  } 
{  number  of  checkers  open  in  stage  n  > 
of  ‘Intflrrayl;  {  optieal  decision  table  ) 

of  A Rea I Array  1;  {  forward  recursion  values  } 
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L  :  Integer; 
next_L  :  integer; 

R  :  Intflrray2; 
a  :  raol; 
iargat  :  raat; 
teopl  :  integer; 
taap2  :  raol; 
taap3  :  raol; 

R  :  raol; 

ij  :  intagar; 

dap  :  intagar; 

tiaa  :  array IDoyOfUeeki  of 

CurrentDay  :  OayOfUaak; 

FirstDay  :  DoyOfUaak; 

LastOay  :  OayOfUaak; 

choica  :  intagar;  <  usad 

dona  :  boo loan; 

FirstDate  :  data; 

LastOata  ;  data; 
CurrantOata  :  data; 
txcBBQfiOiaS  poo  I  ®on  | 


{  nuabar  of  customers  in  I ina  } 
{  nuabar  of  customers  in  I ina  at  naxt  staga  > 
{  nuabar  of  arrivals  in  stags  n  } 
{  custoaar  aaiting  tiaa  } 
{  das i rad  custoaar  aaiting  tiaa  } 
{  taaporory  variable  > 
{  taaporory  variable  > 
{  taaporory  variable  > 
{  return  function  value  } 
{  intagar  loop  control  variables  } 
{  nuabar  of  custoaar  departures  in  stage  n  > 
hours;  {  dai ly  hours  of  the  store  > 

{  currant  day  of  the  week  } 
{  1st  day  of  scheduling  period  } 
(  last  day  of  scheduling  period  } 
in  eenu  to  select  portion  of  prograe  to  run  > 
{  terminates  the  program  } 
{  1st  date  of  scheduling  period  > 
{  last  date  of  scheduling  period  ) 
{  currant  data  in  scheduling  period  } 
{  tails  if  MaxHours  mas  exceeded  } 


<« 


function  calc-nun^cust  <last_jxaLjcust,arrivals,nua_sarvers  ; 

integer;  interval ,serv ice-rate  :  real)  :  integer; 

< 

This  function  calculates  the  nuabar  of  customers  in  line  at 
the  end  of  the  next  stage  given: 

last-nuB-cust  ■  •  customers  at  the  beginning  of  the  stage 
arrivals  ■  •  customers  arriving  during  the  stage 
num-servors  *  *  servers  open  during  the  stage 
service-rate  *  mean  service  rate  of  all  open  servers 
interval  *  length  of  time  in  the  staga 

> 

var 

teap  :  integer;  {  store  *  customers  to  check  if  it  is  nonnegative  > 
begin 

teap  :*  last-nuaucust  +  arrivals  -  trunc< interval *sarvica_rate%*jm-servers> 
if  teap  <  0  then 
calc-niaucust  :*  0 

else 

ealc-num^cust  :■  teap; 

end. 


function  ma i t i ng_t i me  <numjcust,next_numjcust,num_servers  : 

integer;  service. time  :  real)  :  real; 

( 

This  function  calculates  the  aeon  customer  malting  time  in  time 
period  n  given: 

nueucust  *  number  of  customers  in  line  at  the  beginning  of 
time  period  n 

next_nuB_jcust  ■  number  of  customers  in  I  irte  at  the  end  of 
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tiM  period  n 

nus-servers  ■  nusber  of  servers  open  during  tine  period  n 
serv ice-tiee  ■  mcti  service  tiee  of  oil  servers 

> 

begin 

saiting-tiee  :■  <nus_cust  +  next_nus_cust)  *  0.9  *  service-ties  /  nus-servers 
end; 

function  Return  <Uait,desiredjeait  :  real)  :  real; 

{ 

This  function  is  the  return  function  of  the  dynasic  progressing 
forsulation: 

Rlnl  «  I  Uaitlnl  -  desired_jsai  t  I 

•here: 

Uait  ■  the  eean  custoser  salting  tise 
desiredjeait  *  desired  custoser  salting  tise 

> 

begin 

Return  :*  RbsCUait  -  desiredjeait); 
end; 

{sseeessweeeeeees  m  m  mmmmmmm  >m  n  nim  n  i  \t*mm*m****<mmm* } 

procedure  Deallocate; 

{ 

If  the  dynaefc  progrossing  algor i the  results  in  an  optisai 
allocation  of  checkers  shich  is  above  the  lisit  on  total 
hours,  this  procedure  is  used  to  resove  checkers  froe  stages 
so  that  the  change  in  the  return  function  is  Mini si zed. 

) 

var 

stage  :  integer; 

best-place  :  integer;  <  best  stage  to  resove  a  server  > 

z  :  real;  {  current  best  value  of  objective  function  } 

nes-z  :  real;  {  value  of  objective  function  to  be  tested  } 

begin 


shite  (TotalHours  >  MaxHours)  do 
begin 

z  :»  IHFIHITV; 

far  n  :■  I  to  HN  do  {  find  best  place  to  resove  a  server  > 

begin 

if  kCnl  >  HINjCHECKERS  then 
begin 

k(nl  :■  kin]  -  I;  {  reduce  nusber  of  servers  at  stage  n  } 

{  Calculate  the  nee  return  } 

CurrentDay  :■  FirstDay; 
hour  :■  tise (CurrentDay] .open; 
nes_z  :■  0; 


for  stag*  :*  1  to  Nrl  do 
begin 

if  (hour  ■  (tioo l CurrontDay]. close  -  int(TN/80.0*t00>>>  than 
next_L  :«  0 

•Iso 

rwxLL  :•  calc_nuojcust<L,fl[stago],klstago],TN,MU>; 

•  :»  oaiting_tiM<L,noxt-L,k[stago],MUINU); 

R  :»  Re  tum<e,  TARGET); 
nooj:  :■  noe_x  ♦  R; 

If  (hour  ■  (tioo I CurrontDay]. close  -  int(TN/60.0*100>>>  thon 
begin 

if  CurrontDay  ■  Sunday  than 
CurrontDay  :*  Monday 

•Iso 

CurrontDay  :*  succ(CurrontOay); 
hour  :»  tioel CurrontDay], open; 

«nd 

•iso 

hour  :■  hour  +  int(TM/00. 0*100); 
ond;  {  for  stago  :■  1  to  NN  > 


{  !f  tho  noo  rotum  is  better,  roooobor  it  } 

if  noo-z  <  z  thon 
bog  in 

z  :«  noo_z; 
bos t_p loco  :■  n; 
end;  {  if  noo_z  <  z  } 

ktnl  :■  klnl  ♦  1;  {  rostoro  nuobor  of  servers  at  stago  n  } 

•nd;  {  if  klnl  >  MIMXHECKERS  > 

•nd;  {  for  n  :■  1  to  HM  } 

{  Rooovo  server  at  tho  boot  placo  } 
ktbost-plocol  :«  ktbost_ploco)  -  1; 
ond;  {  o*>i lo  (Total Hours  >  MaxHours)  > 
ond;  {  procedure  Doal locate  } 


procedure  Fi i IDPTablo; 
oar 

ch  :  char; 


begin 

CloarScroon; 

•ritoInC Starting  to  fill  OP  state-stage  table* >; 
CurrontDay  :»  LastDay; 

if  ti •• [CurrontDay ] . open  ■  tioo [CurrontDay). close  thon 


iff 


begin 

if  CurrentOay  ■  Monday  then 
CurrentOay  :*  Sunday 

else 

CurrentOay  :■  pred(CurrontDay); 

end; 


hoir  :»  tiee (CurrentOay] .close;  {  Begin  with  last  stage  ) 

n  :«  MM; 

eritelnCn  -  *,HM>; 
for  L  :■  0  to  LJIflX  do 
begin 

teep2  :■  100000.0; 

if  «TM  *  MU  *  MRXjCHECKERS)  <  L)  then 
begin 

teepl  :■  HflX_CHECKERS: 

e?»  eai ting_tiee<L,0,MRX-CHECKERS,MUINU); 

R  :■  Retum(e, target); 
teep2  :■  R; 

end  {  if  «TN  *  MU  *  MRXJCHECKERS)  <  L>  } 

else 

begin 

for  j  :■  M IN-CHECKERS  to  MRXJCHECKERS  do 
begin 

if  <<TN  *  MU  ♦  j)  >»  <L  +  RJMMJ))  then 
begin 

•  :■  eaiting_tiee(L,OJ,MUIMU>; 

R  :■  Retumte,  target); 
if  R  <  teep2  then 
begin 

teep2  :■  R; 
teepl  :*  j; 
end;  {  if  R  <  teep2  } 
end;  {  if  <<TN  *  MU  *  j )  >-  CL+RIMM1 >>  ) 
end;  {  J  :«  MIMJCHECKERS  to  MRXJCHECKERS  } 
end;  {  else  > 
d(NH]*[U  :*  teepl; 
f (NN1* (LI  :»  teap2; 
end;  <  for  L  :»  0  to  LMRX  } 
hour  :■  hour  -  lnt(TN/5Q.0*100>; 

for  n  :■  MM-1  down to  1  do 
begin 

•ritelnC'n  «  ’,n>; 

if  (hour  ■  ti  eel  CurrentOay].  open)  then 
begin 

teep2  100000.0; 

for  J  n IN-CHECKERS  to  MRX-CHECKERS  do 

begin 

next-L  :■  calc_rK«ucust<0,flln],  j  ,TN,MU); 
if  next-L  <  LMRX  then 
begin 

m  :■  eai  ting-tiec<0, next-L ,j  ,MUIMU>; 

R  :»  Re  tum<»,  target); 
if  <R+f ln+1 J* Inext-Ll )  <  teep2  then 
begin 

teep2  :■  R+f(n+H*lnext-LJ; 


I 
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tMpt  j; 

• 

•nd;  {  If  iNxti  <  LHRX  th «n  } 

•nd;  {  j  HINjCHECKERS  to  HRXjCHECKERS  > 
din)* 10)  :■  taepl; 
f Ini* 101  :«  taap2; 

•nd  <  if  hour  *  tieetCurrentDayJ.open  } 

•  lM 

begin 

if  (hour  ■  (tieelCurrantDagl.close))  then 
begin 

for  L  :»  0  to  LMRX  do 
begin 

teep2  :»  100000.0; 

if  «TN  *  HU  *  HRXjCHECKERS)  <  <L  +  fitnl»  then 
begin 

teepl  :■  MflXJCHECKERS; 

next-L  :■  calc-jiuejiXBtCL^ftln^MRXjacCKERS^^nU) 
•  «ai t ing_ti«*<L ,r*xt_L, HRXjCHECKERS, MU INU>; 

R  :*  Retumte,  target); 
teep2  R  ♦  ftneiriOl; 

end  {  if  «TN  *  HU  *  HRXjCHECKERS >  <  <L  ♦  Rlnl»  } 

else 

begin 

for  j  :•  H IN-CHECKERS  to  HRXJCHECMERS  do 
begin 

if  «TH  ♦  HU  *  j)  >■  <L  +  fllnl»  then 
begin 

•  :■  eaiting_tiee(L,0,j,HUINU>; 

R  :■  Retum<», target); 
if  «R  ♦  fln+tl*101)  <  teep2)  then 
begin 

teep2  :■  R  ♦  fln+ll‘10); 
teepl  :■  j; 

•nd;  {  if  R  ♦  f In+1 J* 10)  <  t«ep2  } 
end;  {  if  «TN  *  HU  *  j)  >■  L)  } 
end;  {  j  :«  HINJCHECKERS  to  HRXJCHECKERS  } 
end;  {  else  } 
dtnriU  :«  teepl; 
ftnl'lU  t«ep2; 
end;  {  for  L  :■  0  to  LMRX  ) 

•nd 

else  <  normal  hour  of  the  dag  } 
begin 

for  L  :■  0  to  LHflX  do 
begin 

t«ep2  :■  100000.0; 

for  j  :■  HINjCHECKERS  to  HRXJCHECKERS  do 
begin 

next_L  :«  calc-Jxeucu*t<L,RInJJ,TN,HU>; 
if  nextj.  <  LHflX  then 
begin 

•  :»  eaiting_tiee<L,next-LJ,HUINU>; 

R  :■  Retum<e, target); 
if  <R+ftn+1J*(next_Ll>  <  teep2  then 
begin 

teep2  :■  R+ftn+D‘lnextJ.1; 
teepl  :■  j; 
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and;  {  If  (fi+f (n+irinaxt_Ll>  <  tasp2  > 
and;  {  if  naxt_L  <  LMHX  > 
and;  (  j  (UN-CHECKERS  to  HRXJCHECKERS  > 
dlnl'lL!  :■  taapl; 
flnlMLl  :■  ta*>2; 
and;  {  L  0  to  LMRX  ) 
and;  {  also  > 
and;  {  alsa  > 

if  (hour  ■  tiaal Curran tDoyl.opan)  than 
bag  in 

if  CirrantOay  ■  Nondag  than 
Curran (Day  :■  Sunday 

alsa 

CurrantOay  :■  Prad(CurrantDay); 
if  tiaa(CurrantOayl.opan  ■  t iaa [CurrantOay l.c I osa  than 
bag  in 

If  CurrantOay  ■  Nonday  than 
CurrantOay  :■  Sunday 

alsa 

CurrantOay  :■  Prad(CurrantOay); 

and; 

hour  :»  t iaa (CurrantOay l.c I osa; 
and  {  if  (hour  ■  t iaa (CurrantOay l.opan)  > 

alsa 

hour  :■  hour  -  int(TN/60. 0*100); 
and;  {  n  :■  NN-1  down to  t  } 
and;  {  procadura  FillDPTabla  } 

procadura  ForaardPass; 
var 

ch  :  chcr; 
bag  in 

«ritaln( 'Starting  forward  pass  in  DP'); 

L  0; 

CurrantOay  :*  FirstOay; 

if  tiae (CurrantOay l.opan  *  tiaal Curran tDay 1 .closa  than 
bagin 

i f  CurrantOay  *  Sunday  than 
CurrantOay  :■  Monday 

alsa 

CurrantOay  :■  succ(CurrantOay); 

and; 

hour  :■  t iaa (CurrantOay 1 .opan; 

Total Hours  :»  0; 

for  n  :■  1  to  NN  do 
bagin 

aritalnCn  ■  '  ,n>; 
k(nl  :■  d(nl*(Ll; 

TotalHours  :»  TotalHours  ♦  k(nl; 
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if  (hour  ■  (tiaatCirrantDay 1  closa))  than 
noxt_L  :»  0 

•Im 

naxt-L  :■  calc-Jiua_cust(L,fl[n],k{n],TN,MU); 

L  :■  naxt-L; 

if  (hour  ■  tiaa  (CurrantOay].  closa)  than 
bag  in 

if  (CurrantOay  ■  Sunday)  than 
CurrantOay  :■  Monday 

alsa 

CurrantOay  :■  Succ(CurrantOay); 
if  tiaa [Curran tOay].opan  ■  tiaa (CurrantOay]. cl 
bag  in 

if  CurrantOay  ■  Sunday  than 
CurrantOay  :■  Monday 

alsa 

CurrantOay  :•  succ (CurrantOay); 


hour  :■  tiaa l CurrantOay] .opan; 
and  {  if  (hour  *  t i aa ( Curran tDayl.c I osa)  ) 

alsa 

hour  :■  hour  ♦  int(TN/CO.  0*100); 
and;  {  n  :■  1  to  MM  } 
ana, t  procaoura  roraararass  / 


procodura  Gatflrri vals(var  MuaStagas: intagar;  war  arrivals: lntArray2>; 


war 

InputFila  :  taxt; 
n  :  intagar; 

bagin 

rasatdnputFi  la,  'flrrivals.dat' ); 
n  :■  0; 

•hi  la  not  SaakEofdnputFi  la)  do 
bagin 

n  :*  n  ♦  1; 

if  not  SoakEoIndnputFi  la)  than 
raoddnputFi  la,arriwals(nl) 

alsa 

bagin 

raodln; 

raoddnputFi  la,  arrival  sin]) 
and; 

and; 

ClosadnputFf  la); 

NuaS togas  :■  n; 


procadura  ImtHours; 

{ 

This  procadura  sats  dafoult  stora  hours  as: 


Dog 

Open 

Close 

Hondag 

0000 

0000 

Tuesday 

0900 

2000 

Wednesday 

0900 

1900 

Thursday 

0900 

2100 

Friday 

0900 

1900 

Saturday 

0800 

1800 

Sunday 

0900 

1700 

war 

dag  :  DayOfUeek; 
bag  in 

ties (Monday!. open  :»  0000; 
tiee (Monday!. close  :«  0000; 
ties [Tuesday], open  :■  0000; 
tiee (Tuesday!. close  :«  2000; 
ties (Wednesday!. open  0900; 
t i m (Wednesday 1 . c I ose  :■  1000; 
tiee (Thursday 1. open  :*  0000; 
tiea(Thursdayl. close  :■  2100; 
tiee (Friday!. open  ;■  0000; 
tiae(Fridayl. close  :■  1000; 
tiee (Saturday!. open  :■  0800; 
t  iM  (Saturday  1.  cl  as*  .*1800; 
tiee(Suidayl.open  :•  0000; 
tiee (Sunday!. cl ose  :■  1700- 
end;  {  procedure  InitHours  } 


procedure  Ini  tPointers<nueber_jof_stages  :  integer); 

This  procedure  initializes  the  pointers  for  dynaeic  eeeory 
al location. 

> 

begin 

for  n  :»  1  to  nueber_of_s togas  do 
begin 
Nee<d(nl>; 
hee<f(nl); 
end; 

end; 

procedure  InitTarget; 

{ 

This  procedure  sets  the  default  custoeer  waiting  tiee  goat 
to  5.0  einutes . 

) 

begin 

target  :■  5.0; 

end;  {  procedure  InitTarget  } 


procedure  InltfkwHours; 

This  procedure  sets  the  default  eaxieue  olloeable  checker  hours 
to  infinity. 

} 

bog  in 

haxHours  :•  INTJMFINITY; 
end;  {  procedure  Ini thaxMours  } 


(♦eeeeeemi  in  h  i  hi  m»»hh  »ii<  ii  1 1  m  eeeeeeeee *****  m  i  »»»♦»»»»►»»♦»» 

procedure  henu <var  selection  :  integer); 

begin 

Cl ear Screen; 

GotoXV<23, 10  >; 

•riteCI.  Run  Phase  I*); 

0otoXV<23, I1>; 
erite('2.  Run  Phase  If); 

QotoXY(23, 12); 

*rite<'3.  Run  Phase  I  t  Phase  ll‘>; 

0otoXY<23, 13 >; 

erite<’4.  Set  Weekly  Store  Operating  Hours); 

0otoXV<23, 14); 

erite<*3.  Set  Desired  Custoeer  Waiting  Ties'); 

GotoXY<23, 13); 

erite<’6.  Set  Lieit  on  Total  Checker  Hours'); 

GotoXY(23, 10); 
erite<‘7.  Quit); 

0otoXV<23, 19); 

eri te< 'Enter  Selection:*); 

read I n<se I ect i on ); 

Cl ear Screen; 
end; 


nm  i  imiiii  weeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 

procedure  ReadDate<var  FirstOate,  LastOate  :  date); 

This  procedure  asks  the  user  for  the  first  and  last  day  to 
be  scheduled  and  converts  the  input  into  the  proper  date  format. 


oar 

FirstDay,  LastOay  :  string!  101; 

BadResponse  :  boolean; 

begin 

Cl ear Screen; 
repeat 

BadResponse  :*  false; 

GotoXYC 12, 12); 

eri te< ‘Enter  first  day  of  scheduling  period  (ee/dd/yyyy): ' ); 
readln<FirstOay); 

if  <ord(First0ayI11)  >  49)  or  <ord<First0ay[1))  <  48)  then 
BadResponse  :■  true; 


if  <ord<First0ayt21)  >  37)  or  <ord<FirstDay[21>  <  48)  then 
BadResponse  :*  true; 


if  (FirstOayM)  ■  *  V >  and  <ord<First0ayC21)  >  SO)  than 
BadRasponsa  :■  trm; 

if  <ord<FirstDay(41)  >  51)  or  <ord<FlrstDay[41>  <  48)  than 
BojOipflm  ;»  trua* 

if  <ord<FirstDay(51)  >  57)  or  <ord<FirstDay[51>  <  48)  than 
BadRasponsa  :■  trua; 

if  (ord<FirstOagI7))  >  50)  or  <ord<FirstDaym)  <  40)  than 
BadRasponsa  ;*  trua; 

if  <ord<FirstDay(8I )  >  57)  or  <ord<FlrstDayl8]>  <  48)  than 
BadRasponsa  :*  trua; 

if  <ord<First0ay{01)  >  57)  or  <ord<FirstDayt91>  <  48)  than 
BadRasponsa  :■  trua; 

if  (ard(First0ag(101)  >  57)  or  <ord(First0ay(101)  <  48)  than 
BodRasponsa  ;*  trua' 

if  <First0ayI3i  <>  •/>  or  <FirstDayt6)  <>  '/')  itm\ 

BadRasponsa  :■  trua; 

if  BadRasponsa  than 
bagin 

ClaarScraan; 

OotoXV< 12, 10); 

aritaf Incorrect  format  or  out  of  possible  range,  try  again.'); 
and 


all 


bagin 

GotoXV<  12, 10); 
ClaarEOL;  • 


until  not  BadRasponsa; 
rapaat 

Dodnasponsa  ; *  falsa; 

GotoXYC 12, 13); 

arita< ‘Enter  last  day  of  scheduling  period  (aa/dd/yyyy): * >; 
raadln<LastOay); 

if  Cord<LastDoy[1])  >  49)  or  (ord(LastOayM))  <  48)  than 
BadRasponsa  :■  trua; 

57)  or  (ord<LastOayI21)  <  48)  than 


if  <ord<LastPayI21)  > 
BadRasponsa  :•  trua 
if  <LastBay( 1]  «  '1') 
BadRasponsa  :*  trua 
if  <ord<LastOay (4 ] )  > 

DQOnapOnn  .*  uMI 
if  <ord(Last0ayl31 >  > 
BadRasponsa  :■  trua 
if  (ord<Last0ay(71)  > 
BadRasponsa  :■  trua 
if  <ord<Last0ay(81)  > 
BadRasponsa  ;■  trua 
if  <ord<Last0ay(9]>  > 
BadRasponsa  :■  trua 
if  <ord<LastOay(IO)) 
BadRasponsa  :■  trua 
if  <Last0ayC31  <>  '/' 
BadRasponsa  :»  trua 
if  BodRasponsa  than 
begin 

QotoXV<61, 13); 

ClaarEOL; 

OotoXV< 12,10); 


and  <ord(Last0ayI2)>  >  50)  than 
St)  or  <ord<LastOay(41)  <  48)  than 
57)  or  <ord(LastOay(5})  <  48)  than 
50)  or  <ord(LastOay[?I)  <  48)  than 
57)  or  <ord(LastDay(8])  <  48)  than 
57)  or  <ord(Last0ay(9])  <  48)  than 
>  57)  or  < ord ( Last Day [ 10J)  <  48)  than 
)  or  <LastDayt&]  <>  '/')  than 
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mriteC Incorrect  format  or  out  of  possible  range,  try  again. ’ >; 
end 

else 

begin 

GotoXV<  12,10); 

CtmarCOL; 

•nd; 

until  not  BadResponse; 

{  Convert  date  froe  string  format  to  'date*  format  (integer)  > 

FtrstOate. month  :•  10  *  (ord(FirstDaglll)  -  48)  ♦  <ord<FlrstOayI21)  -  48); 
Firs tOate. dag  :■  10  *  (ord(First0ag[41 )  -  48)  +  (ord(FirstOagtSl)  -  48); 
FirstOate.gear  :■  1000  *  (ord(First0ag[71>  -  48)  ♦  100  *  <ord<FirstDay(8))  - 
48) 

♦  10  *  (ord(FirstDagl9) )  -  48)  ♦  (ord(FirstOagtlOl)  -  48); 

LastOate. month  :■  10  *  (ord(LastOagMl)  -  48)  ♦  <ord<LastDay(2J)  -  48); 
LastOate.dag  :■  10  *  (ord(Last0ag(41)  -  48)  ♦  (ord<LastOag(31 >  -  48); 
LastOate. gear  :■  1000  *  <ord<Last0ag(7I)  -  48)  ♦  100  *  <ord<Last0ag(81 )  -  48) 

♦  10  *  (ord(LastOag(91)  -  48)  ♦  (ordCLastDagMOl)  -  48); 

end;  {  procedure  ReadDate  > 


function  Jul i an< SomeOate  :  date)  :  integer; 

This  function  takes  a  date  in  'standard'  date  format  and 
concerts  it  to  Julian  format  (  gggg  and  ddd  ) 


war  . 

temp  :  integer; 
begin 

case  SomeOate. month  of 


1 

temp 

■  SomeOate. dag; 

2 

temp 

■  SomeOate. dog  *  31; 

3 

temp 

*  SomeOate. dag  +  39; 

4 

temp 

■  SomeOate .dag  ♦  90; 

3 

temp 

*  SomeOate .dag  ♦  120, 

6 

temp 

■  SomeOate .dag  *  151 

7 

temp 

-  SomeOate. dog  +  181 

8 

temp 

■  SomeOate. dag  ♦  212 

9 

temp 

■  SomeOate. dag  ♦  243 

10 

temp 

■  SomeOate. dag  ♦  273 

11 

temp 

■  SomeOate. dag  +  304 

12 

temp 

■  SomeOate. dog  ♦  334 

end;  {  case  SomeOate .month  of  } 
if  ((SomeOate. gear  mod  4)  ■  0)  then 


temp  :■  temp  ♦  1; 
Julian  :*  temp; 
end;  {  function  Julian  > 


function  0ags8etmeen(First0ate, LastOate  :  date)  :  integer; 


This  function  calculates  tha  nuabar  of  days  bote— n  too  dates. 

}  i 

war 

i  intagor; 

FirstJUl  ian,  LasUul  ian  :  integer; 
toop  :  intogor; 

bog in  | 

FirstJul ian  :■  Jul ian<FfrstOato); 

LastJulian  :■  Jul ian(LastOata); 
toop  :■  LastJulian  -  FirstJul ian; 
for  I  :■  FirstOato.yoar  to  LastOate.yoar-1  do 
if  <1  ood  4)  -  0  thon 

toop  :■  toop  +  366  ■ 

«lM  I 

toop  :■  toop  +  363; 

OaysOeteoon  :■  toop; 
end;  {  function  OaysSotooon  > 


procedure  SetDateCwar  StartOay,EndOay  :  DayOfUook;  war  NuoStag—  :  integer); 

This  procedure  asks  the  user  far  the  first  and  last  day  to 
bo  scheduled.  It  teen  calculates  tee  total  nuobor  of  stages 
(hours)  in  the  schedul ing  period  and  the  day  of  the  seek  for  the 
first  and  last  days. 


war 

NumDay*,temp  :  integer; 

Ref  :  date; 

OayRrray  :  array [0. .61  of  OayOfUeek; 

today  :  OayOfUeek;  | 

begin 

Ref.eonth  :■  1;  {  Reference  date  is  Friday,  1  January  1988  } 

Ref. day  :■  1; 

Ref. year  :■  1988; 

OayRrraylOl  :*  Friday;  I 

OayRrrayMl  :«  Saturday;  1 

0ayRrray(21  :*  Sunday; 

OayRrrayO]  :*  Monday; 

0ayRrray(41  :■  Tuesday; 

Dayf¥ray(3]  :*  Wednesday; 

Doyflrrayl6J  :■  Thursday; 

RoadDato(F I rs tOa to, Las tOa to ); 

(  Determine  the  day  of  the  meek  for  the  first  day  to  bo  scheduled  } 

NuoOays  :■  Days8eteeen(Ref,FirstDate); 
toop  :■  NuoOays  ood  7; 

StartOay  :■  OayRrray I toop I ;  I 

{  Determine  the  day  of  the  meek  for  the  last  day  to  be  scheduled  } 

NuoOays  : *  DaysBoteeen  <Ref,LastDate); 
toop  :■  NuoOays  ood  7; 
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EndDoy  :*  DoyArrayltaapl; 


NueDays  :■  Days8etween<FfrstOate,LastOate); 
today  :»  StcrtOay; 

NuoStagos  :■  0; 
for  teop  :»  0  to  NuoOays  do 
bog  in 

if  <t loot today l.opon  <>  tiseC today] .close)  than 
NuoStages  :»  NuoStages  ♦  t loot today! .close  -  t loo (today l.opon  ♦ 
int<TN/G0. 0*100); 

if  today  *  Sunday  than 
today  :■  Monday 
olso 

today  :»  Succ(  today); 

•nd; 

NuiStQQts  :»  NusSiagM  div  100; 
if  TN  <>  GO  than 

NuoStages  NuoStages  *  round<60/TN>; 
and;  {  procoduro  SotOato  } 


procoduro  MlnCostFlow<NueNodes,Nueflrcs  :  intogor; 

c  :  flrcflrray; 
war  x  .  nr  array, 

F  :  flrcflrray; 

T  :  flrcflrray; 

b:  NodoArrau: 
u  :  flrcflrray); 

This  procoduro  solves  the  Miniouo  Cost  Floo  probloo: 


Min  C‘x 
s.t.  fix  ■  b 
0  i  x  S  u 


k 


x  *  vector  of  arc  flows 
c  *  vector  of  arc  floo  costs 
b  *  vector  of  node  supply  or  desands 
u  ■  vector  of  arc  floo  upper  bounds 
A  *  node  arc  incidence  oatrix 

NuoNodes  *  nueber  of  nodes  in  the  network  (sax  is  254) 

NuoArcs  ■  nueber  of  arcs  in  the  network  (wax  is  255  -  MueNodes) 

Notice  that  the  node  arc  Incidence  oatrix  was  not  one  of  the 
input  paraoeters .  Instead,  the  arrays  F  t  T  ere  used,  where: 

F  ■  “Froo"  function  of  an  ere,  i.e.  FtxC I , j )1  ■  i 
T  ■  “To"  function  of  an  arc,  i.e.  T(x<i,j)1  »  j 
Notice  that  FIT  require  ouch  less  eeeory  than  A. 

The  following  constants  oust  be  defined  before  the  procedure 
is  cal  led: 

Maxflrcs  *  255; 

MaxNodes  *  254; 

The  following  types  oust  be  defined  before  the  procedure 
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it  cal  ltd: 

flrcflrray  *  array M.  .haxflrcs  1  of  intagar; 

Nodaflrray  ■  array! 1. .ftaxHodas)  of  intagar; 

For  addad  «p— d,  ail  variables  and  para— tart  ara  dafinad  at 
intagart.  If  nonintagar  values  ara  needed,  tiaply  radafina 
tha  variables  and  paraMtart  at  raalt. 

Tha  out  of  ki I  tar  algori tha  it  utad  to  solve  tha  problaa. 
Rafaranca:  Rigor i that  for  Hataork  Progressing,  by  Jeff  L. 

(Canning tor  t  Richard  U.  Haigaton.  John  Uilay 
t  Sort,  Now  York,  I960. 


const 

inf  -  9999; 
typa 

SlackNodaSat  ■  tat  of  1..naxNodas; 
SlackflrcSat  ■  tat  of  1..Maxflrcs; 


war 

pi  :  NodaArray; 

cost  :  Rrcflrray; 

thata  :  intagar; 

ptil  :  tat  of  1..f1axArcs; 

pti2  :  sat  of  1..naxRrcs; 

Nhat  :  tat  of  1. .Haxttadat; 

flhat  :  sat  of  1..Haxflrcs; 

dalta  :  Nodaflrray; 
i  :  intagar; 
j  :  intagar; 
s  :  intagar; 

naa<Lxrc  :  boo loan;  {  indicatas 


{  dual  variables  > 
{  cost  l  j  1  «  pilFIjll  -  pi  IT  I  j  11  -  cljl  } 
{  aoount  to  change  dual  variables  } 
{  candidate  sat  for  traa  } 
{  candidate  sat  for  traa  } 
{  currant  nodas  in  traa  > 
{  currant  arc*  in  traa  } 
{  oat  to  change  float  in  cycle  } 
{  loop  control  variable  > 
{  loop  control  variable  } 
{  out  of  ki Iter  arc  > 
ohathar  an  out  of  ki I  ter  arc  is  knoan  } 


in-kiltar  :  boolean;  {  indicates  ohathar  all  arcs  ara  in  kilter  } 
no-cycla  :  boolean;  {  indicatas  ohathar  a  cycle  exists  } 
NusSlackRrcs  :  intagar;  {  total  •  of  arcs  after  slack  arcs  ara  added  } 
NuaS I ackhodas  :  intagar;  {  total  •  of  nodas  aftar  slack  node  is  addad  } 


procadura  Ini  tiaUSolution; 

{ 

This  procedure  finds  an  initial  feasible  solution  to  start 
tha  algor i the. 


Tha  initial  solution  is  found  by  use  of  tha  “al l-arti f icial 
start.*  fin  extra  node,  huahodas+1,  ft  addad  to  tha  network. 
For  each  source  node  i  with  supply  blil,  a  slack  arc  is 
addad  frow  i  to  NusModas+1  with  cltlock  arc)  ■  0, 
ultlack  arc)  ■  infinity,  and  xtslack  arc)  «  b(i). 

For  each  dawand  node  k  with  dasand  |b(kl|,  a  slack  arc  is 
addad  from  MusNodas+1  to  k  with  ctslack  arc)  *  infinity, 
u [slack  arc)  *  infinity,  and  xlslack  arc)  »  Iblkl). 


} 


ML  :  Mt  of  1.  .NoxNodes; 

MU  :  Mt  of  1.  .NoxNodes; 
found  :  boo loan; 
j  :  integer; 
i  :  integer; 

bog  in 

for  j  :■  1  to  Nuoflrcs  do 
bog  in 
xtj 1  :■  0; 
orvd; 

for  i  :»  1  to  NuoNodes  do 
bogin 

if  bC i I  >0  thon 
bogin 

xti+Nuoflrcsl  :■  bti  1; 
ct i+Nuoflrcsl  :■  0; 
uti+Nueflrcsl  :■  inf; 

Tli+Nueflrcsl  :■  NuoModos  ♦  1; 

Ft i+Nueflrcsl  :»  i; 
ond;  <  if  bill  >  0  } 
if  bti 1  <*  0  thon 
bogin 

xti+Nueflrcsl  :•  -blil; 
cti+Nuoflrcsl  :■  inf; 
uti+Mueflrcsl  :■  inf; 

Tti+Muoflrcsl  :■  i; 

Fli+Mueflrcsl  :■  NuoModos  +  !; 
ond;  {  if  bill  <»  0  } 
ond;  {  for  i  :■  1  to  NuoModos  } 

for  i  :»  1  to  NuoSlackNodos  do  {  fill  pi  til  for  ail  nodas  > 

pitil  0; 

for  j  :•  1  to  NuoSlockRrcs  do  {  cooputa  costtjl  for  all  arcs  } 
cost  l  j  1  :«  pilFIJll  -  pi  IT  I  j  1 1  -  ctj  ); 

ond;  {  procoduro  Initial-Solution  } 


I 


procoduro  flrcSearch<var  i n_k i I  tor  :  boolean;  var  os  :  integer );  - 

< 

This  procoduro  soarchs  for  an  out  of  kilter  arc  in  the  network. 

If  all  network  arcs  are  in  kilter,  thon  the  flag  in-kilter  is 
Mt  to  true.  Otherwise,  the  first  out  of  kilter  arc  found  is 

returned  in  es,  and  the  flag  in-kilter  is  Mt  to  faiM.  1 

This  procoduro  assuees  that  all  pgr awe tors  are  available,  ie 
ctjl,  xCJ  1,  pi  (nl,  t  costtjl  ■  pi  IF  I  j  11  -  pitTIjll  -  ctjl 

bogin 

in— ki I  ter  :»  true;  . 

j  0;  5 

repeat 

j  :■  J  ♦  I; 

if  (costtjl  <  0)  and  <xtjl  >  0>  thon 
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InJUItar  falsa; 

if  (costtjl  >  0)  and  (xljl  <  ut j 1 >  than 
in-kilter  :■  falsa; 

until  <not  iruKiltar)  or  <j  »  MueSlackflrcs); 
if  (not  in_kilter>  than 


procedure  prieal (as  :  integer;  war  no-cy-le-f  lag  :  boolean); 

This  procedure  executes  the  priaal  phase  of  the  out 
of  kilter  algorithe  eith  arc  as.  If  the  algor i the 
terei nates  with  the  conclusion  that  no  cycle  exists, 
r*o_cyc I e_f I ag  is  set  to  true.  0 there i se, 
no_cyc I e_f I og  is  set  to  false. 

} 


war 

found  :  boolean; 
j,i  :  integer; 
tree-path  :  flrcflrray; 
nue_pathjarcs  :  integer; 
current-node  :  integer; 

procedure  f ind_path(k, I  :  integer;  M  :  SlackhodeSet; 
fl  :  SlackflrcSet; 
war  path  :  fircflrray; 
war  nue_ in-path  :  integer); 

This  procedure  finds  the  path  thru  the  tree  N 
from  node  k  to  node  I  (assuees  that  there  is 
only  a  single  path  in  N  froe  k  to  I).  The  path 
is  returned  in  the  array  path  and  the  nueber  of 
arcs  in  the  path  is  returned  in  nue_in_path. 

> 

war 

next-node  :  array ( 1. .ttaxhodesl  of  Integer; 
found-nodes  :  SlackModeSet; 
i,j,  current-node  :  integer; 
foirtd  :  boolean; 

begin 

next-nodell!  :*  I;  {  start  at  I  -  there’s  no  next  node  } 

found-nodes  :■  I  I  I; 

repeat 

for  J  :■  1  to  NueSlockflrcs  do 
begin 

if  (j  in  fl)  then 
begin 

if  (TIj)  in  found-nodes)  and 
<not<F(jl  in  found-nodes))  then 
begin 

found-nodes  :■  found-nodes  ♦  (Fiji  ); 
next_nodeIFIj ]]  TIj  1; 


•nd; 

if  (Fiji  in  foundLnodas)  and 

<not<T(j]  in  found-nodas>)  than 
bag  in 

found-nodas  :■  found-nodas  +  I  Ttj 1  1; 
noxt-nodatTljll  FIJI; 

#nd  ’ 

and;  {  if  <Tlj  1  in  N>  and  (Fiji  in  N>  > 
and;  {  for  J  :«  1  to  NuaSlackflrcs  > 
until  (k  in  found_nodas>; 
nua_in_path  :■  0; 
currant_noda  :■  k; 
rapaat 
j  :■  0; 

found  :*  falsa; 
ahila  not  found  do 
bag  in 

j  :■  j  ♦  1; 

if  <<j  in  fl  >  and  <FIj 1  ■  naxt_nodaIcurrant_nodal) 
and  <T I j  I  ■  curran t_noda))  than 
bag  in 

nua_irupath  :»  nua_i n_pa th  +  1; 
found  :»  trua; 
pathlruua-irupathl  :«  j; 
currant_noda  :■  FIJI; 
and; 

if  «j  in  fl>  and  (Fiji  •  currant_noda>  and 
<TIj 1  »  nax t-nodal curran t-nodal ))  than 
bagin 

nua_in_path  :*  nua_in_path  +  1; 
found  :■  trua; 
pathlnua_in_pathl  :*  j; 
curranUnoda  :■  TtjJ; 
and; 

and;  {  ahi la  not  found  } 
until  currant_noda  *  I; 
and;  {  procadura  find_path  } 

{*»a*a*aaaaa*a**aa**aaaaa*a*a*»a*a»aaaaaaa*aaaaaa»aa*a»*a*aaaaaaa*n »♦*♦*♦} 
bagin 

Rhat  :■  II;  {0.  Initializa  } 

if  costtasl  <  0  than 
bagin 

Nhat  :■  I  Tlasl  1; 
daltalTIasll  :■  xlasl; 
and  {  if  cost Iasi  <  0  } 
alsa  {  cost ias]  >«  0  } 
bagin 

Nhat  :■  l  Ftesl  1; 
daltalFIasIl  ulas]  -  xlasl; 
and;  {  alsa  > 

rapaat 

psi  1  :■  11;  {1.  Oataraina  Candidatas  for  Traa  } 

psi2  :*  II; 

f or  J  :■  1  to  NuaSlockflrcs  do 
bagin 


A 

i 


a 


4 


J 


J 


J 
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if  <<j<>es)  and  (costtjl  >■  0>  and  <xtj]  <  ul j  1  >  and 
<TIj  I  in  Hhat>  and  <not<Flj]  in  Nhat)))  than 
psi 1  :■  psi 1  ♦  (  j  1; 

if  <<J<>es)  and  <costIj]  <■  0)  and  (xljl  >  0)  and 
<Flj  3  in  Nhat)  and  <not<Ttj]  in  Nhat)))  than 
psi2  psi2  ♦  I  j  1; 
and;  {  for  j  :«  t  to  NuaSlackArcs  } 
if  <<psi 1  +  psl2)  •  (1>  than 
bag  in 

no_cycle_f lag  :■  trua; 
exit; 

and  {  if  ((psil  +  psi2)  «  U>  > 

else 

no_cycla_f lag  :*  falsa; 

j  ;*  0;  {2.  Append  Noe  Arc  to  Tree  } 

found  :■  falsa; 
ropaot 
j  j  ♦  t; 
if  (j  in  psil)  than 
bagin 

if  (daltalTIj 13  <  <uf j 1  -  xlj 1 >>  than 
daltalFIj ]]  :■  daltalTIj]] 

else 

daltalFljll  :■  (uljl  -  xljl); 
found  :•  trua; 
and; 

if  (j  in  psi2>  than 
bagin 

if  (daltalFIj ]]  <  xljl)  then 
daltalTIj]]  daltalFIj]] 

else 

daltalTIj  1]  :*  xljl; 
found  :»  true; 
end; 

until  found; 

Nhat  :*  Nhat  +  I  Tlj |,FIJ 1  1; 

Ahat  :■  Ahat  +  I  j  1; 

until  <FIes]  in  Hhat)  and  <Tlesl  in  Nhat); 

if  (costlesl  <  0)  then  {  3.  Breakthrough  } 

begin 

f  ind_path<TIesl,Fles),Mhat,Ahat,  tree-path, nun-path-arcs); 
currant-node  :*  Ties]; 
for  j  :*  I  to  nua-path-orcs  do 
begin 

if  <T l tree-path l j 11  ■  current-node)  then 
begin 

xttree-pathlj 1)  :»  xltree-pathlj ]]  ♦  deltalFlesll; 
current-node  :*  Fltree_pathlj  11; 
end 

else  {  FI  tree-pa thtj ]]  3  current-node  > 
begin 

xltree-pathlj]]  :*  xltree-pathlj IJ  -  del talFIes] ]; 
current-node  :■  Tttree-pathlj  II; 
end; 

end;  {  for  j  :*  1  to  nua_path_arcs  } 
xfssl  :■  xles]  -  deltalFlesll; 


end 

else  {  costlesl  >■  0  } 


begin 

current-node  :■  Fie*]; 

f  ind_path<F(esl,Tlesl,Nhat,Ahat,  tree-path, nue-path-jarcs); 
for  j  :■  1  to  nue_path_jarcs  do 
begin 

if  <T I tree-path I j 11  ■  current-node)  then 
begin 

xl  tree-pa  thlj  11  :■  x(  tree-pa  thlj  1]  +  deitatFIe*]]; 
current-node  :■  Fttree-pathlj 11; 


else  {  FI  tree-path £j 11  *  current-node  } 
i  d 

xltree-pathlj ]]  :*  xl tree-pa thlj II  -  deltalFIesll; 
current-node  :■  T I  tree-pa thlj 13; 
end; 

end;  {  for  j  :*  1  to  nue-path_rrcs  > 
x(es)  xles)  ♦  del ta (Ties 3 3 ; 
end; 

end;  {  procedure  prieal  > 


procedure  dual (es  :  integer;  uar  need-arc-f lag  :  boolean); 

{ 

This  procedure  iepleeents  the  dual  phase  of  the  out  of 
kilter  algorithm  on  the  tree  developed  in  the  prieal 
phase  and  sets  the  flag  need-arc-f  lag.  If  arc  es  is  in 
kilter  at  the  end  of  the  dual  phase,  a  nee  out  of  kilter 
arc  oust  be  found,  so  need-arc-f lag  is  set  to  true.  If 
arc  es  is  still  out  of  kilter  at  the  end  of  the  dual  phase, 
then  the  prieal  phase  oust  be  executed  again,  and  the 
need-arc-flag  is  set  to  false. 

> 


var 

i,j  :  integer; 
begin 

{  Start  with  T  ■  (Mhat.flhat)  >  {  0.  Initialization  > 

{  developed  in  prieal  phase  } 

psi 1  :»  13;  {  1.  Detereine  Arcs  Incident  on  T  } 

psi2  :*  11; 

for  j  :■  1  to  MueSlackflrcs  do 
begin 

If  (cost  I  j  3  <  0)  and  (noUFIjl  in  Mhat>>  and 
(TIj  3  in  Nhat)  then 
psi 1  :■  psi  1  +  l  j  1; 
if  (cost! j  1  >  0)  and  (FIJI  In  Mhat)  and 
<not(TIj  3  in  Mhat)>  then 
psi2  :■  psi2  +  I  j  3; 
end;  <  for  j  :■  1  to  NuoSlackflrcs  } 

theta  :»  inf;  {  2.  Detereine  haxieue  Pereissable  Change  } 

for  j  :■  1  to  NueSlackArcs  do 
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if  (j  in  (psil  +  psi2>>  and  (Abs(cost[j 1)  <  theta)  than 
theta  :■  flbsCcostlj 1); 


for  i  :>  1  to  NumSIackNodes  do  {  3.  Reduce  Duals  } 

if  <  i  in  Mhat  >  then 
pi C i  ]  pill!  -  theta; 

for  j  :*  1  to  MueSlackflrcs  do  {  Update  Costs  > 

costfjl  :■  pitFIjll  -  pi  ITI j  1]  -  ctj  1; 

need  nrc_f lag  :«  true;  {  Set  needier  c-f lag  ) 

if  < cos ties]  <  0)  and  (xlesl  >  0)  then 
niid  nr r  firm  :>  fain* 
if  (cos ties]  >  0)  and  (xlesl  <  u(esl)  then 
neecLjarc-flag  :■  false; 

end;  {  procedure  dual  > 


begin 

NumSIackNodas  NueNodes  ♦  1; 

MueSlackflrcs  :«  Numflrcs  +  huehodes; 

Initial-Solution;  {  Find  an  initial  set  of  feasible  floes  } 

need-arc  :■  true;  {  Do  not  currently  knoe  an  out  of  kilter  arc  } 

repeat 

if  need-arc  then  {  Find  an  out  of  kilter  arc  s  } 

begin 

flrcSearch<in-ki lter,s>;  {  If  no  out  of  kilter  arcs, 

in-kilter  ■  true  } 

end; 

if.  not  in-kilter  then 


begin 

primal <s,no_cycle);  {  Execute  the  prieal  phase  eith  arc  s  } 

need  arc  :■  true; 
if  <no-cycle)  then 
begin 

dual <s, need-arc);  {  If  prieal  phase  finds  no  cycles, 

execute  dual  phase  } 

end; 

end;  {  if  not  i n_k i I  ter  } 
until  in_kilter; 
end;  {  procedure  hinCostFloe  } 


procedure  Phasel; 

{ 

This  procedure  per fores  the  Phase  I  portion  of  the  problem,  i.e. 
determination  of  the  checker  requirements. 

The  checker  requirements  are  returned  in  the  global  vector  k  and 
are  output  to  the  file  Servers. out. 

} 

var 

OutputFile  :  text; 


begirt 


S*tDat*<Fir»tDay,Le»tDog, Total Stages);  {  Asks  for  period  to  schedule  } 
GetArrivals(MN,R);  {  Read  arrival  data  into  A,  *  stages  into  HN  > 


if  (NN  <>  TotaiStages)  then  {  if  NN  <>  TotalStages,  soee thing  wrong  > 
begin 

ClearScreen; 

OotoXV< 13, 12); 

erite(‘Nuaber  of  stages  in  Arrivals.Dat  does  not  Batch  number'); 
GotoXV(15, 13); 

eriteCof  stages  froe  first  to  last  day  of  scheduling  period.’ >; 
OotoXV< 13, 13); 

erite<‘  Strike  return  to  continue. ' ); 

GotoXV<23, 17); 
eriteCNN  •  ’,Hh); 

GotoXV<25, 18); 

eri te( ‘TotalStages  ■  ', Total Stages); 
readln; 

Exit; 


InitPointers(NH);  {  Initial  ize  pointers  } 

FillOPTabfe;  {  Calculates  f  t  d  for  the  entire  OP  state-stage  table  ) 
ForwardPass;  {  Finds  optieai  path  thru  OP  state-stage  table  > 

if  TotalHours  >  NaxHours  then 
begin 

ExceedHours  :■  true; 

Deal  locate; 
end 

else 

ExceedHours  :*  false; 

reerite<OutputFi I e, ‘Servers. out' );  {  Open  output  file  > 

case  FirstOay  of  {  Write  FirstDay  to  output  file  > 

Monday  :  eri teln(OutputFf le, ‘Monday’ ); 

Tuesday  :  writeln<OutputFi le, 'Tuesday* ); 

Wednesday  :  sri teln<OutputFi le, ‘Wednesday ); 

Thursday  :  eri tel n((kitputFi ie, 'Thursday* ); 

Friday  :  eri ta!n<OutputFi le, ‘Friday’ >; 

Saturday  :  eri teln(OutputFi le, ‘Saturday’ ); 

Sunday  :  eri tel n<OutputFi I e, 'Sunday' ); 
end;  {  case  FirstOay  } 

for  n  :*  1  to  NN  do  {  Write  optieai  •  checkers  to  ‘Servers. out’  } 

eri teln(OutputFi le,k(nl ); 

CloseCOutputFi le);  {  Close  output  file  } 

end;  {  procedure  Phase I  } 


procedure  Phase  1 1 ; 

This  procedure  reads  the  first  day  of  the  scheduling  period 
and  the  ideal  checker  requirements  (output  from  Phase  I) 
froe  the  file  ' Servers. out' .  It  then  determines  the  optimal 
number  of  checkers  to  schedule  for  each  shift.  These  optimal 
shifts  are  output  to  the  file  ‘Shifts. out' . 
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InputFi le  :  text; 

OutputFile  :  text; 

OayText  :  string! 10); 

MuaHours  :  integer; 

MuaShifts  :  integer; 

NueVars  :  intagar; 

Nua6  :  intagar; 

Nun?  :  intagar; 

HuaO  :  intagar; 

Hurt  :  intagar; 

Hua4  :  intagar; 
shifUF  :  flrcflrroy; 
shi ft_T  :  flrcflrroy; 
shift  :  flrcflrroy; 
b  :  Modeflrray; 
obj  flrcflrroy; 
uppar  :  flrcflrroy; 
shift_start  :  intagar; 
shift~and  :  intagar; 

bag  in 

resetdnputFl  I e, 'Servers. out'  );  {  Opan  Input  Fila  ) 

readlndnputFi le, OayText);  {  Rood  first  day  of  scheduling  period  } 
if  OayText  ■  'Monday'  then 
FirstOay  :■  Monday; 
if  OayText  *  'Tuesday'  then 
FirstOay  :«  Tuesday; 
if  OayText  ■  'Wednesday'  then 
FirstOay  :*  Uedrtesday; 
if  OayText  ■  ‘Thursday*  then 
FirstOay  :■  Thursday; 
if  OayText  ■  ‘Friday*  than 
FirstOay  :■  Friday; 
if  OayText  ■  ‘Saturday*  then 
FirstOay  :■  Saturday; 
if  DayText  *  ‘Sunday*  then 
FirstOay  :■  Sunday; 


{  File  variable  for  'Servers. out’  > 
{  Fite  variable  for  ‘Shifts. out'  } 

(  •  of  hours  store  is  open  for  current  day  > 
{  *  possible  daily  shifts  for  current  day  } 
{  *  possible  daily  shifts  +  *  deviation  vars  > 
{  *  possible  8  hour  shifts  for  current  day  ) 

<  •  possible  7  hour  shifts  for  current  day  > 

{  •  possible  0  hour  shifts  for  current  day  } 

{  •  possible  3  hour  shifts  for  current  day  } 

{  *  possible  4  hour  shifts  for  current  day  } 

{  “Froe*  functions  for  arcs  (shifts)  } 
{  “To"  functions  for  ares  (shifts)  > 
{  *  servers  in  each  shift  } 
{  change  in  hourly  requireeents  } 
{  cost  coefficients  in  objective  function  } 
{  upper  bounds  on  shift  variables  } 
{  start  tiae  for  the  current  shift  } 
{  end  tiae  for  the  currant  shift  } 


{  Read  checker  requireeents  into  k  and  total  •  stages  into  HH  ) 
n  :■  0; 

ehile  not  SeekEof( InputFi le)  do 
begin 

n  :«  n  ♦  1; 

if  not  SeekEo I n( InputFi le)  then 
read( InputFi I  a, kin]) 

else 

begin 
read In; 

readdnputFi  le,ktnl) 


end; 

ClosednputFi  le); 

MM  :»  n; 

retri te(OutputFi le, 'Shifts. out' >; 
if  ExceedHours  then 
begin 

«ri te(OutputF) le, ' **e*aa**ee* 


••); 
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eritaln<OutputFI  le,  '***  iwwwwmmxmw^miwwh'  >; 

mri  taCOutputFi  le,  'URRHINO  URRHINO  URRNINQ  '); 

mri teln(OutputFi  le,  ‘URRNINQ  URRNINQ  URRNINQ  URRNINQ' >; 

mritaln(OutputFile); 

mri  telnCOutputFi  le); 

mrite<OutputFi le, '  Ideal  checker  requirements  for  each  hour  ' >; 
eriteln<OutputFi le, ‘exceed  the  eaxieue  total  hours.'); 
mriteln(OutputFi le); 

erlte<OutputFI le, '  Suggest  you  reduce  target  customer  malting  ' ); 
mri teln(OutputFi le, ‘time  and  run  program  again.'); 
mriteln<OutputFi le); 
mriteln<OutputFi le); 

mriteCOutputFi le, 'URRNINQ  URRNINQ  URRHINO  '); 

sri teln<OutputFi le, 'URRHINO  URRHINO  URRHINO  URRHINO'); 

mri teln<Outputf i le); 

eri  te<OutputFi  le, 1 »  mmmmmmmmmmmmmmmm  ■  >; 

mri  te(n<OutputFi  le,  * ****mm*mmm*mm mm *«* i ****** ******* •  >; 
mri  teln<OutputFi  le); 
mriteln<OutputFi  le); 
end; 


CurrentOay  :■  FirstOay; 
n  :»  1; 


repeat 


NumHours  :■  (time! CurrentOay]. close  -  time (CurrentOay] .open)  div  100; 
if  NumHours  <>  0  then 
begin 

NumHours  :■  NumHours  +  1;  {  Add  1  hoir  for  additional, 

redundant  node  } 


for  I  :■  1  to  NumHours  do 
begin 

if  i  »  1  then 
b(i )  :■  kin) 

else 

if  i  ■  NumHours  then 
bill  :>  -klnl 

else 

bill  :*  kin)  -  kln-1); 
if  i  <  NumHours  then 
n  :■  n  ♦  1; 


Num8  NumHours  -  8; 

Nun?  :■  NumHours  -  7; 

Num6  :■  NumHours  -  0; 

Num5  :■  NumHours  -  9; 

Num4  :■  NumHours  -  4; 

NumUars  :■  NumO  *  Num?  +  Num6  +  Num5  +  Num4  ♦ 
<2  *  (NumHours  -  1)); 

for  j  1  to  NumB  do 
begin 

shift-FIJl  j; 

shift-Ttj ]  j  ♦  8; 

objljl  0; 
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upper! j 1  HRXJCHECKERS; 

end; 

for  j  :«  (Hurt  +  1)  to  (Hurt  *  Hue7>  do 
begin 

shift-Fljl  :■  j  -  Hurt; 
shlft_TlJ )  J  -  Hue6  ♦  7; 
obj  IJ  1  0; 

upper (J I  :•  HRXJCHECKERS; 

* 

for  j' :■  (Hurt  +  Hue?  +  1)  to  (Hurt  +  Hue?  +  Nue6>  do 
begin 

*hi ft— Ftj 1  j  -  (Hue©  +  Hue?); 

shift-Tljl  :■  j  -  (Hue©  ♦  Hue?)  ♦  6; 
obj  Ij  J  0; 

upper (j ]  :*  HRXJCHECKERS; 

for  j':*  (Hue©  ♦  Hue?  +  Hurt  ♦  1)  to 

<Hue8  ♦  Hue?  +  Hurt  +  Hurt)  do 


begin 

shi ft F [ j  I  :«  j  -  (Hue©  +  Hue?  ♦  Hurt); 

shlft_Ttj 1  :■  j  -  (Hurt  +  Hue?  ♦  Hurt)  ♦  3; 
obj  Ij )  :«  0; 

upper ij 1  .-HRXJCHECKERS; 
end; 

for  j  (Hurt  ♦  Hue?  ♦  Hurt  ♦  Hurt  ♦  1 )  to 

(Hurt  ♦  Hue?  +  Hurt  ♦  Hurt  +  Hue4)  do 


begin 

shift-Fljl  :«  J  -  (Hurt  ♦  Hue?  +  Hurt  +  Hurt); 
shift-Tlj 1  :»  j  -  (Hurt  ♦  Hue?  +  Hurt  ♦  Hurt)  +  4; 
objtjl  :«  0; 

upper tj ]  HRXJCHECKERS; 
end; 


{  “Froe"  t  "To"  functions  for  deviation  variables  } 


j  :»  (Hurt  +  Hue?  ♦  Hurt  ♦  Hurt  +  Hue4  ♦  1); 
for  i  :«  1  to  (HueHours-1)  do 


begin 

shift-Fljl  :•  i; 
shift_Ttj  1  i  «•  1; 

obj  IJ  1  I; 

upper (j  ]  :>  HRXJCHECKERS; 
shift-Flj  +  11  :•  i  +  1; 
shift_Tlj+11  :■  i; 
objCj+11  :■  1; 

(jppertj  +  1 1  HRXJCHECKERS; 

j  :■  j  +  2; 

end; 


{  di-  } 

{  di+  } 

< 


HinCostFloe(HueHours,MueUars,obj,shift,shift-F,shift_T,b,upper>; 


{  Urite  shifts  to  output  file  > 
wri teln(0utputFi le); 

case  CurrentDag  of  j 

Honday  :  eriteln(OutputFi le, ’Monday’ >; 

Tuesday  :  srlteln(OutputFi le, ’Tuesday' >; 

Wednesday  :  srl te I n(0utputFi I e, ’Wednesday’ >; 

Thursday  :  eri teln(0utputFi le, ‘Thursday’ ); 
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Friday  :  oritalntOutputFi la, 'Friday'  >; 

Saturday  :  oritalntOutputFi  la, 'Saturday' ); 

Sunday  :  oritalntOutputFila, 'Sunday' >; 
and;  {  com  CurrantOay  of  } 

casa  CurrantOay  of 
Monday  :  ari ta I n< 'Monday' >; 

Tu— day  :  or  I  taint' Tuesday '  ); 

Wednesday  :  ori taint 'Wednesday' ); 

Thursday  :  or I  taint 'Thursday' >; 

Friday  :  ori taint 'Friday' >; 

Saturday  :  ori taint 'Saturday '); 

Sunday  :  ori taint 'Sunday' >; 
and;  (  casa  CurrantOay  of  > 

oritalntOutputFi la, '  8  Hour  Shifts: ' >; 

f or  j  :*  1  to  Nua8  do 
begin 

shift-start  :■  tiaalCurrantOayl .open  ♦  ttj  -  1)  *  100); 
shi ft_and  :■  shift_start  ♦  <800), 
ori tatOutputFi la, ‘  *>; 

oritalntOutputFi  la,shi ft-start:4,  ’“',shifl_and:4, 

*  *,shiftlj  1:2); 

and; 

ori talntOutputFi la, '  7  Hour  Shi fts: ' >; 

for  j  :■  <Mua8  +  t)  to  tNuaS  ♦  Hue?)  do 
begin 

shift-start  :■  tiaalCurrantOayl .open  ♦ 
ttj  -  Nua6  -  1>  *  100>; 
shiftjand  :■  shift-start  ♦  <700); 
ori tatOutputFi la, ‘  ' >; 

oritalntOutputFi la, shi ft-start:4, '-’,shift_and:4, 
',shiftlj 1:2); 

and; 

oritalntOutputFi la, '  6  Hour  Shifts: * >; 

for  j  <Mua8  ♦  Mua7  +  1)  to  tNuaS  +  Mua7  +  NuoO)  do 
begin 

shift-start  :*  tiaa I Curran tOayl.opan  + 

ttj  -  Mue8  -  Mua7  -  1)  +  100); 
shift-end  :*  shift-start  +  tOOO); 
or (tatOutputFi I a,'  * >; 

oritalntOutputFi la,shift_start:4, ’-’,shiftjend:4, 

•  ’ ,shi ftlj 1:2); 

and; 

ori talntOutputFi I a,’  5  Hour  Shi fts: ’ ); 

for  j  :*  tNua8  ♦  Mua7  ♦  Nuab  ♦  1)  to 

<Nue8  +  hue?  +  Nua6  +  Nua5)  do 

begin 

shift-start  :»  tiaalCurrantOayl. open  + 

ttj  -  Mue8  -  Mua7  -  Muad  -  1)  *  100); 
shift-snd  :■  shift-start  +  t300>; 
ori tatOutputFi la, 1  ’>; 

ori talntOutputFi la, shift-start: 4, shift-snd: 4, 

' ,shl ftlj 1:2); 

and; 
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■ritelntOutputFlle, 1  4  Hour  Shifts: ' 

for  j  :■  (NueB  +  Nue7  +  MueA  +  Nue5  ♦  1>  to 

(NueB  ♦  Hum?  +  Hurt  ♦  Hurt  +  Mum4)  do 

bog  in 

shift-start  :■  tiee (CurrantOay] .open  ♦ 

«j  -  Hurt  -  Hue?  -  Hurt  -  Hue 5  -  1)  *  100); 
shifljend  :■  shift-start  +  (400); 

■rite(OutoutFi la. *  • >• 

■ritelntOutputFI \ a, shift-start: 4, ,-,#shiftj*nd:4/ 

'  ‘,shift(j  1:2); 

and; 

end;  <  if  HueHours  <>  0  } 

if  CurrantOay  *  Sunday  than 
CurrantOay  :■  Monday 

also 

CurrantOay  :■  succ<CurrantOay); 

unti  I  <n  >■  NN>; 

Closa(OutputFi  la); 
aritaln(>Muabar  of  stages  «  1 ,n>; 
and;  {  procadura  Phase  II  } 

procedure  Solve; 

begin 
Phase I ; 

Phase! I ; 

and;  {  procadura  Solve  } 


procadura  SatHours; 

{ 

This  procedure  is  used  to  Modify  the  store's  hours  of  operation. 


correct,  BadResponsa  :  boo  loan; 
day  :  DayOfUaak; 
rasponsa  :  char; 

begin 

correct  :■  false; 

■hi  la  not  correct  do 
begin 

Cl ear Screen; 

GotoXY<1,7); 

•ritelnC 
■ri tain; 

■riteln; 

■ri taln< ’ 

■riteln; 

for  day  :*  Monday  to  Sunday  do 


{  ClrScr  in  IBM  Turbo  } 
Current  hours  of  operation  are: ' > 

Day  Open  Close'); 


til 


begin 


case  day  of 
Monday 

•riteC 

Monday  ' ); 

Tuesday 

•rlteC 

Tuesday  ’ ); 

Wednesday 

•riteC 

Wednesday' ); 

Thursday 

•riteC 

Thursday  ‘ ); 

Friday 

•riteC 

Friday  ’); 

Saturday 

•rlteC 

Saturday  ' ); 

Sunday 

•riteC 

Sunday  ' ) 

end;  {  case  > 

if  t i •• [day J .open  ■  tiaetdayl .close  then 
•ritelnC  Closed') 

else 

•r  f  teln(  t  i»e(day  ].  open:8,  t  ieelday ).  close:  8>; 
end;  {  dag  *  Mon  to  Sun  } 

•riteln; 

•riteln; 

BodResponse  :•  true; 

•hile  BodResponse  do 
begin 

•riteC  Is  this  correct  (V/H)?’ >; 

read I n<response ); 

if  (response  ■  ‘g'>  or  (response  ■  *V*>  then 
exi  t 

else 

if  (response  ■  'n')  or  (response  ■  'M‘ >  then 
DodTIesponse  :■  false 

else 

eriteln( 'Please  use:  V  for  yes,  N  for  no‘>; 
end;  {  BodResponse  } 

Cl ear Screen; 


•riteln; 

•riteln; 

•ritelnC 
■ritelnC 
•riteln; 

•ritelnC 
•ritelnC 
>; 

■riteln; 

for  day  :*  Monday  to  Sunday  do 
begin 

case  dag  of 
Monday  :  writelnC 
Tuesday  :  •ritelnC 
Wednesday  :  •ritelnC 
Thursday  :  erltelnC 
Friday  :  ■ritelnC 
Saturday  :  •ritelnC 
Sunday  :  ■ritelnC 


Please  use  24-hour  tiees  for  all  entries.’); 
Tiees  eust  be  rounded  to  nearest  hour. 1 ); 


For  days  whan  store  is  closed,  enter  0000' ); 
for  opening  tine  and  0000  for  closing 


Monday* ); 
Tuesday' ); 
Wednesday’ ); 
Thursday ' >; 
Friday’ >; 

Saturday' ); 
Sunday’ >; 


•rlteC 

read  I  n(ti  ee  [day],  open); 
■riteC 

read  in(  t  i  ee  [day  1 .  c  I  ose ); 
end; 

end;  {  not  correct  } 
end;  {  procedure  SetHours  } 


Open  (xxxx):  ’  >; 
Close  (xxxx): ’ ); 


I 

I 

i 
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procedure  SetftaxHours; 

This  procedure  proepts  the  user  for  the  eaxieue  allowable 
checker  hours  and  stores  it  in  the  global  variable  MaxHours . 

begin 

Cl ear Screen; 

OotoXV< 15, 12); 

erite< ‘Enter  eaxieue  allowable  total  checker  hours: '); 
read I nChaxHours ) ; 
end;  {  procedure  SetflaxHours  > 


procedure  SetTarget; 

{ 

This  procedure  proepts  the  user  for  the  desired  cus tower 
waiting  tiee  and  stores  it  in  the  global  variable  target. 

> 

begin 

Cl  ear  Screen; 

GotoXY(20, 12); 

write< ‘Enter  desired  cus tower  waiting  tiee:‘); 
read I n< target); 
end; 


begin 

InitHours;  {  Initialize  the  default  store  hours  } 

InitTarget;  {  Initialize  the  default  desired  custower  waiting  tiee  } 
InitHaxHours;  {  Initialize  the  default  eaxieue  total  checker  hours  } 

done  :*  false; 
repeat 

Henufchoice);  {  Put  up  eenu  } 

case  choice  of 

1  :  Phase); 

2  :  Phased; 

3  :  Solve; 

4  :  SetHours; 

5  :  SetTarget; 

6  :  SetflaxHours; 

7  :  done  :■  true; 
end;  {  case  } 

until  done; 

Exit; 

end. 
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Date 

Day 

Time 

A(n) 

27-May-87 

Wednesday 

9 

240 

27- May-87 

Wednesday 

10 

312 

27- May-87 

Wednesday 

11 

347 

27- May-87 

Wednesday 

12 

329 

27-May-87 

Wednesday 

13 

326 

27- May-87 

Wednesday 

14 

323 

27- May-87 

Wednesday 

15 

319 

27-May-87 

Wednesday 

16 

384 

27- May-87 

Wednesday 

17 

337 

27-May-87 

Wednesday 

18 

131 

28- May-87 

Thursday 

9 

158 

28- May-87 

Thursday 

10 

263 

28- May-87 

Thursday 

11 

311 

28- May-87 

Thursday 

12 

270 

28- May-87 

Thursday 

13 

300 

28- May-87 

Thursday 

14 

292 

28- May-87 

Thursday 

15 

309 

28- May-87 

Thursday 

16 

369 

28- May-87 

Thursday 

17 

337 

28- May-87 

Thursday 

18 

275 

28- May-87 

Thursday 

19 

201 

28- May-87 

Thursday 

20 

57 

29- May-87 

Friday 

9 

62 

29- May-87 

Friday 

10 

153 

29- May-87 

Friday 

11 

170 

29- May-87 

Friday 

12 

227 

29- May-87 

Friday 

13 

258 

29- May-87 

Friday 

14 

291 

29- May-87 

Friday 

15 

316 

29- May-87 

Friday 

16 

322 

29- May-87 

Friday 

17 

278 

29- May-87 

Friday 

18 

114 

30- May-87 

Saturday 

8 

166 

30- May-87 

Saturday 

9 

290 

30- May-87 

Saturday 

10 

345 

30-May-87 

Saturday 

11 

354 

30- May-87 

Saturday 

12 

400 

30- May-87 

Saturday 

13 

350 

30- May-87 

Saturday 

14 

363 

86 


30- May-87 

Saturday 

30- May-87 

Saturday 

30- May-87 

Saturday 

31- May-87 

Sunday 

31-May-87 

Sunday 

31- May-87 

Sunday 

31- May-87 

Sunday 

31 -May-87 

Sunday 

31-May-87 

Sunday 

31-May-87 

Sunday 

31-May-87 

Sunday 

2-Jun-87 

Tuesday 

2-Jun-87 

Tuesday 

2-Jun-87 

Tuesday 

2-Jun-87 

Tuesday 

2-Jun~87 

Tuesday 

2-Jun-87 

Tuesday 

2-Jun-87 

Tuesday 

2-Jun-87 

Tuesday 

2-Jun-87 

Tuesday 

2-Jun-87 

Tuesday 

2-Jun-87 

Tuesday 

3-Jun-87 

Wednesday 

3-Jun-87 

Wednesday 

3-Jun-87 

Wednesday 

3-Jun-87 

Wednesday 

3-Jun-87 

Wednesday 

3-Jun-87 

Wednesday 

3-Jun-87 

Wednesday 

3-Jun-87 

Wednesday 

3-Jun-87 

Wednesday 

3-Jun-87 

Wednesday 

4-Jun-87 

Thursday 

4-Jun-87 

Thursday 

4-Jun-87 

Thursday 

4-Jun-87 

Thursday 

4-Jun-87 

Thursday 

4-Jun-87 

Thursday 

4-Jun-87 

Thursday 

4-Jun-87 

Thursday 

4-Jun-87 

Thursday 

4-Jun-87 

Thursday 

4-J\in-87 

Thursday 

4-Jun-87 

Thursday 

20 

52 

5-Jun-87 

Friday 

9 

157 

5-Jun-87 

Friday 

10 

263 

5-Jun-87 

Friday 

11 

317 

5-Jun-87 

Friday 

12 

286 

5-Jun-87 

Friday 

13 

287 

5-Jun-87 

Friday 

14 

304 

5-Jun-87 

Friday 

15 

319 

5-Jun-87 

Friday 

16 

275 

5-Jun-87 

Friday 

17 

286 

5-Jun-87 

Friday 

18 

96 

6-Jun-87 

Saturday 

8 

120 

6-Jun-87 

Saturday 

9 

235 

6-Jun-87 

Saturday 

10 

281 

6-Jun-87 

Saturday 

11 

345 

6-Jun-87 

Saturday 

12 

340 

6-Jun-87 

Saturday 

13 

334 

6-Jun-87 

Saturday 

14 

338 

6-Jun-87 

Saturday 

15 

308 

6-Jun-87 

Saturday 

16 

315 

6-Jun-87 

Saturday 

17 

129 

7-Jun-87 

Sunday 

9 

1 

7-Jun-87 

Sunday 

10 

244 

7-Jun-87 

Sunday 

11 

295 

7-Jun-87 

Sunday 

12 

280 

7-Jun-87 

Sunday 

13 

312 

7-Jun-87 

Sunday 

14 

353 

7-Jun-87 

Sunday 

15 

333 

7-Jun-87 

Sunday 

16 

124 

9-Jun-87 

Tuesday 

9 

197 

9-Jun-87 

Tuesday 

10 

295 

9-Jun-87 

Tuesday 

11 

303 

9-Jun-87 

Tuesday 

12 

304 

9-Jun-87 

Tuesday 

13 

316 

9-Jun-87 

Tuesday 

14 

278 

9-Jun-87 

Tuesday 

15 

306 

9-Jun-87 

Tuesday 

16 

352 

9-Jun-87 

Tuesday 

17 

309 

9-Jun-87 

Tuesday 

18 

186 

9-Jun-87 

Tuesday 

19 

53 

10-Jun-87 

Wednesday 

9 

119 

10-Jun-87 

Wednesday 

10 

234 

10-Jun-87 

Wednesday 

11 

277 

88 


10-Jun-87 

10-Jun-87 

10-Jun-87 

10-Jun-87 

10-Jun-87 

10-Jun-87 

10- Jun-87 

11- Jun-87 

ll-Jun-87 

ll-Jun-87 

ll-Jun-87 

ll-Jun-87 

ll-Jun-87 

ll-Jun-87 

ll-Jun-87 

ll-Jun-87 

ll-Jun-87 
ll-Jun-87 

11- Jun-87 

12- Jun-87 

12-Jun-87 

12-Jun-87 

12-Jun-87 

12-Jun-87 

12-Jun-87 

12-Jun-87 

12-Jun-87 

12-Jun-87 

12- Jun-87 

13- Jun-87 

13-Jun-87 

13-Jun-87 

13-Jun-87 

13-Jun-87 

13-Jun-87 

13-Jun-87 

13-Jun-87 

13-Jun-87 

13- Jun-87 

14- Jun-87 

14-Jun-87 

14-Jun-87 

14-Jun-87 

► 


Wednesday 

12 

249 

Wednesday 

13 

238 

Wednesday 

14 

275 

Wednesday 

15 

271 

Wednesday 

16 

259 

Wednesday 

17 

246 

Wednesday 

18 

0 

Thursday 

9 

73 

Thursday 

10 

188 

Thursday 

11 

242 

Thursday 

12 

184 

Thursday 

13 

207 

Thursday 

14 

284 

Thursday 

15 

298 

Thursday 

16 

302 

Thursday 

17 

316 

Thursday 

18 

213 

Thursday 

19 

193 

Thursday 

20 

58 

Friday 

9 

173 

Friday 

10 

235 

Friday 

11 

292 

Friday 

12 

276 

Friday 

13 

255 

Friday 

14 

317 

Friday 

15 

307 

Friday 

16 

301 

Friday 

17 

310 

Friday 

18 

97 

Saturday 

8 

75 

Saturday 

9 

189 

Saturday 

10 

313 

Saturday 

11 

328 

Saturday 

12 

363 

Saturday 

13 

361 

Saturday 

14 

327 

Saturday 

15 

366 

Saturday 

16 

355 

Saturday 

17 

145 

Sunday 

9 

68 

Sunday 

10 

236 

Sunday 

11 

299 

Sunday 

12 

339 

89 


14-Jun-87 

Sunday 

13 

326 

14-Jun-87 

Sunday 

14 

336 

14-Jun-87 

Sunday 

15 

343 

14-Jun-87 

Sunday 

16 

159 

16-Jun-87 

Tuesday 

9 

200 

16-Jun-87 

Tuesday 

10 

322 

16-Jun-87 

Tuesday 

11 

300 

16-Jun-87 

Tuesday 

12 

327 

16-Jun-87 

Tuesday 

13 

279 

16-Jun-87 

Tuesday 

14 

302 

16-Jun-87 

Tuesday 

15 

371 

16-Jun-87 

Tuesday 

16 

334 

16-Jun-87 

Tuesday 

17 

313 

16-Jun-87 

Tuesday 

18 

199 

16-Jun-87 

Tuesday 

19 

59 
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1  GEN,CAPT  FREV,0  6  JUN  87, 11/16/88, t,V,N,V/V,H,V/S,72; 

2  LIMITS,  1,4,500; 

3  NETWORK 

4  RESOURCE /CHECKERS(Q ), 1 ; 

5  CREATE, XX< 1  >, ,  1; 

8  RCTIUITV, ,TNOU.LE. 600,  OK; 

7  RCTIUITV,, TNOU.GT. 600, KIL; 

8  KIL  TERM; 

0  OK  ASSIGN,ATRIB<2>«5. 156; 

10  ASSIGN,RTRIB<3>-ATRIB<1>+ATRIB<2>; 

11  ASSIGN,  RTRIB<4Mf«J<1>; 

12  AUA I T<1), CHECKERS/ 1; 

13  RCT I U I TV/ 1 , RTR I B<2 >j  SERUICE  TIME 

14  FREE, CHECKERS/ 1; 

15  EUENT,  1; 

16  COLCT, INT< 1 ),TIME  IN  SVSTEI1; 

17  COLCT,  INTO), TIME  IN  QUEUE; 

TERM; _ 

CREATE, 60, 60,, 11; 

EUENT, 2; 

RCTIUITV,, TNOU.GE. 60  .AND.  TNOU.LT. 120, P2; 
RCTIUITV, , TNOU.GE. 120  .AND.  TNOU.LT. 190, P3; 
RCTIUITV, , TNOU.GE. 180  .AND.  TNOU.LT. 240, P4; 
RCTIUITV,, TNOU.GE. 240  RNO.  TNOU.LT. 300, P5; 
RCTIUITV,, TNOU.GE. 300  .RNO.  TNOU.LT. 360, P6; 
RCTIUITV,, TNOU.GE. 360  .RNO.  TNOU.LT. 420, P7; 
RCTIUITV,, TNOU.GE. 420  .RNO.  TNOU.LT. 480, P8; 
RCTIUITV,, TNOU.GE. 480  AND.  TNOU.LT. 540, P9, 
RCTIUITV,, TNOU.GE. 540  .ATC.  TNOU.LT. 600, P 10; 
RCTIUITV,, TNOU.GE. 600  .RNO.  TNOU.LT. 660, P1 1; 
ACT  I U I TV, , TNOU . GE . 660, P 12 ; 

P2  RSSIGN,XX<1 >-600/235.0; 

ALTER, CHECKERS/* 11; 

TERMINATE; 

P3  ASSIGN, XX < 1 >>60.0/281.0; 

ALTER, CHECKERS /+3, 

TERMINATE, 

P4  RSSIGN,XX<1 >60.0/345.0; 

ALTER, CHECKERS/*?; 

TERMINATE; 

P5  RSSIGN,XX<1 >60.0/340.0; 

ALTER, CHECKERS/- 1; 

TERMINATE; 

PC  RSSIGN,XX<1 >-60.0/334.0; 

TERMINATE; 

P7  ASSIGN,XX( 1 >>60.0/338.0; 

ALTER, CHECKERS/*  1; 

TERMINATE, 

PS  ASSIGN, XX<1 >-60.0/306.0; 

ALTER, CHECKERS/-5; 

52  TERMINATE; 
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P9  ASS  I  ON,  XX(1  >60. 0/3 15.0; 

ALTER, CHECKERS /+3; 

TERM I MATE; 

P10  RS8I0M,XX<1 >60.0/129.0; 

ALTER, CHECKERS/- 15; 

TERMINATE; 

P11  ALTER, CHECKERS/- 12; 

TERMINATE; 

P12  ASS  I  ON,  XX ( 1  >100000; 

TERMINATE; 

ENONETUORK; 

INTLC,XX<  1  >0 . 5,  XX<2  >0,  XX<3  >0,  XX<4  >0,  XX<5>0,  XX<6  >0 ; 
INIT,0.0,800.0,N; 

FIN; 
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SLAM  I  I  SUMMARY  REPORT 


SIMULATION  PROJECT  D  6  JUN  87  BV  CRPT  FREV 

DATE  11/16/1088  RUN  NUMBER  1  OF  1 

CURRENT  TINE  0.8000E+03 

STATISTICAL  ARRAYS  CLEARED  AT  TINE  0  OOOOE+OO 

♦•STATISTICS  FOR  VARIABLES  BASED  ON  OBSERVATION** 

MEAN  STANDARD  COEFF.  OF  MINIMUM  MAXIMUM  NO. OF 

VALUE  DEVIATION  VARIATION  VALUE  VALUE  OBS 

TIME  IN  SYSTEM  0.787E+01  0.11IE+01  0.141E+00  0.516E+01  0. 124E+02  2743 

TIME  IN  QUEUE  0.271E+01  0.111E+01  0.409E+00  O.OOOE+OO  0.722E+01  2743 

**FILE  STATISTICS** 

FILE  AVERAGE  STANOARO  MAXIMUM  CURRENT  AVERAGE 

NUMBER  LABEL/TYPE  LENGTH  DEVIATION  LENGTH  LENGTH  WRIT  TIME 

1  AWAIT  0.304  7.508  28  0  2.714 

2  CALENDAR  19.504  12.193  33  1  0.931 

♦♦REGULAR  ACTIVITY  STATISTICS** 

ACTIVITY  AVERAGE  STANDARD  MAXIMUM  CURRENT  ENTITY 

INOEX/LABEL  UTILIZATION  DEVIATION  UTIL  UTIL  COUNT 

1  SERVICE  TIME  17.0787  11.9310  30  0  2743 

♦♦RESOURCE  STATISTICS** 

RESOURCE  RESOURCE  CURRENT  AVERAGE  STANOARO  MAXIMUM  CURRENT 

NUMBER  LABEL  CAPACITY  UTIL  DEVIATION  UTIL  UTIL 

1  CHECKERS  1  17.68  11.931  30  0 

RESOURCE  RESOURCE  CURRENT  AVERAGE  MINIMUM  MAXIMUM 

NUMBER  LABEL  AVAILABLE  AVAILABLE  AVAILABLE  AVAILABLE 

1  CHECKERS  1  0.2713  -15  9 
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SIMULflTlon  PROJECT  Con  lamda  1 


BV  CRPT  FREY 


DRTE  11/16/1988 


RUM  HUMBER  50  OF  50 


CURRENT  TIME  0.8000E+03 

STRTISTICRL  RRflflYS  CL ERRED  RT  TIME  O.OOOOE+OO 


♦♦STATISTICS  FOR  VARIABLES  BRSED  OM  OBSERVATION** 

MERM  STRMORRO  COEFF.  OF  MINIMUM  MRXIMUM  HO. OF 

URLUE  OEU I RT I OM  URRIRTIOM  URLUE  URLUE  OBS 

TIME  IM  SYSTEM  0. 117E+02  0. 102E+02  0  074E+OO  0. I55E+01  0.230E+03  **** 

TIME  IM  QUEUE  0.655E+01  0. 100E+02  0. 153E+01  O.OOOE+OO  0.227E+03  **** 


•♦FILE  STATISTICS** 


FILE 

AUERAGE 

STRMORRO 

MAXIMUM 

CURRENT 

AVERAGE 

HUMBER 

LABEL /TYPE 

LENGTH 

OEU 1 RT 1 OM 

LENGTH 

LENGTH 

WAIT  TIME 

1 

AWAIT 

12.801 

12.971 

52 

0 

3.764 

2 

CALENDAR 

19.155 

11.770 

33 

1 

0.889 

♦♦REGULAR  ACTIVITY  STATISTICS** 

ACTIVITY 

AVERAGE 

STRMORRO 

MAXIMUM 

CURRENT 

ENTITY 

INDEX /LABEL 

UTILIZATION 

DEVIATION 

UTIL 

UTIL 

COUNT 

1  SERVICE  TIME 

17.6153 

11.6713 

30 

0 

2721 

♦•RESOURCE  STATISTICS** 

RESOURCE  RESOURCE  CURRENT  AUERAGE  STRMORRO  MRXIMUM  CURRENT 

NUMBER  LABEL  CAPACITY  UTIL  OEU I RT I  OH  UTIL  UTIL 

1  CHECKERS  1  17.62  11.671  30  0 


RESOURCE  RESOURCE  CURRENT  AVERAGE  MINIMUM  MRXIMUM 

HUMBER  LABEL  AUAILABLE  AVAILABLE  AVAILABLE  AVAILABLE 


Case  3:  6  Jun  87 

SLAM  II  sunitflflv  REPORT 

SIMULATION  PROJECT  Con  Lcnda  1  BV  CAPT  FREV 

DATE  11/16/1988  RUN  NUMBER  50  OF  50 


CURRENT  TINE  0.8000E+03 

STATISTICAL  AAAAVS  CLEARED  AT  TINE  0  OOOOE+OO 


♦•STATISTICS  FOR  UARIABLES  BASED  ON  OBSERUAT I ON** 

MEAN  STANOARO  COEFF.  OF  MINIMUM  NRXINUN  NO. OF 

UALUE  DEU I  AT  I  ON  UARIATION  UAL UE  UALUE  DBS 

TINE  IN  SYSTEM  0. 131E+02  0. 123E+02  0.940E+00  0. 159E+01  0.260E+03  **** 

TINE  IN  QUEUE  0 . 798E+0 1  0 . 122E+02  0 . 153E+0 1  0 . OOOE+OO  0 . 256E+03  **** 


♦♦FILE  STATISTICS** 

FILE  AUERAGE  STANDARD  MAXIMUM  CURRENT  AUERAGE 

NUMBER  LABEL/TYPE  LENGTH  DEU I AT I ON  LENGTH  LENGTH  UAIT  TIME 


1 

AUAIT 

20.763 

21.055 

85 

0 

6.014 

2 

CALENORR 

19.519 

12.029 

33 

1 

0.894 

♦♦REGULAR  ACT  1 U 1  TV  STATISTICS** 

ACTIVITY 

AUERAGE 

STANDARD 

MAXIMUM 

CURRENT 

ENTITY 

INDEX /LABEL 

UTILIZATION 

DEVIATION 

UTIL 

UTIL 

COUNT 

1  SERVICE  TIME 

17.6009 

11.6214 

30 

0 

2762 

♦•RESOURCE  STATISTICS** 


RESOURCE 

NUMBER 

RESOURCE 

LABEL 

CURRENT 

CAPACITY 

AVERAGE 

UTIL 

STANOARO  MAXIMUM 
DEVIATION  UTIL 

CURRENT 

UTIL 

1 

CHECKERS 

1 

17.60 

11.621  30 

0 

RESOURCE 

RESOURCE 

CURRENT 

AVERAGE 

MINIMUM  MAXIMUM 

NUMBER  LABEL  AURILABLE  AUAILABLE  AUAILABLE  AUAILABLE 
1  CHECKERS  1  1.2555  -15  18 
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Case  4;  8  Jun  87 

slam  ti  summary  report 

S I MULflT I OM  PROJECT  Con  Lcada  1  BV  CRPT  FREY 

ORTE  11/16/1988  RUM  NUMBER  50  OF  50 


CURRENT  TIME  0.8000E+03 

STRTISTICRL  RRRRVS  CLERRED  RT  TIME  O.OOOOE+OO 


••STATISTICS  FOR  URRIRBLES  BASED  OM  OBSERURT I OM** 

MERM  STRMORRO  COEFF.  OF  MINIMUM  MAXIMUM  MO. OF 

VALUE  OEUIRTIOM  URRIRTIOM  URLUE  VALUE  OBS 

TIME  IN  SYSTEM  0. 172E+02  0. 178E+02  0. 103E+01  0. 155E+01  0.285E+03  **** 

TIME  IN  QUEUE  0. 121E+02  0. 177E+02  0. 147E+01  O.OOOE+OO  0.282E+03  **** 


♦•FILE  STATISTICS** 


FILE 

RUERAGE 

STRMDARO 

MAXIMUM 

CURRENT 

RUERAGE 

NUMBER 

LABEL /TYPE 

LENGTH 

OEUIRTIOM 

LENGTH 

LENGTH 

URIT  TIME 

1 

flURIT 

42.949 

25.762 

92 

7 

12.554 

2 

CALENOAA 

19.729 

11.764 

33 

2 

0.914 

••REGULAR  RCTIUITY  STATISTICS** 

RCTIUITY 

RUERAGE 

STANDARD 

MAXIMUM 

CURRENT 

ENTITY 

INOEX /LABEL 

UTILIZRTIOM 

OEU 1  AT  1 ON 

UTIL 

UTIL 

COUNT 

1  SERVICE  TIME 

17.5882 

11.5355 

30 

1 

2729 

••RESOURCE  STATISTICS** 


RESOURCE 

NUMBER 

RESOURCE 

LABEL 

CURRENT 

CAPACITY 

AVERAGE  STRNDARO  MAXIMUM  CURRENT 
UTIL  DEVIATION  UTIL  UTIL 

1 

CHECKERS 

1 

17.59  11.535  30 

1 

RESOURCE 

NUMBER 

RESOURCE 

LABEL 

CURRENT 

AVAILABLE 

AVERAGE  MINIMUM 

AVAILABLE  AVAILABLE 

MAXIMUM 

AVAILABLE 

1 

CHECKERS 

0 

0.9292  -15 

21 
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SIMULATION  PROJECT  0  11  JUN  87 


BV  CRPT  FREV 


ORTE  11/16/1988 


RUN  NUMBER  1  OF  1 


CURRENT  TINE  0.8000E+03 

STATISTICAL  ARRAYS  CLEARED  AT  TIME  O.OOOOE+OO 


•♦STATISTICS  FOR  UARIA6LES  BASED  ON  OBSERVATION** 

MEAN  STAfCARO  COEFF.  OF  MINIMUM  MAXIMUM  NO. OF 

UALUE  DEVIATION  VARIATION  VALUE  VALUE  OBS 

TIME  IN  SYSTEM  0.876E+01  0. 175E+01  0. 199E+00  0.516E+01  0.167E+02  2558 

TIME  IM  QUEUE  0.361E+01  0.  175E+01  0.484E+00  O.OOOE+OO  0.115E+02  2558 


••FILE  STATISTICS** 


FILE  AVERAGE  STANOARO  MAXIMUM  CURRENT  AtCRAGE 

NUMBER  LABEL /TYPE  LENGTH  DEVIATION  LENGTH  LENGTH  WAIT  TIME 


1 

2 


AWAIT  11.535  7.857  33 

CALENDAR  18.386  8.817  30 


0  3.608 

1  0.955 


••REGULAR  ACTIVITY  STATISTICS** 

ACTIVITY  AVERAGE  STANDARD  MAXIMUM  CURRENT  ENTITY 

INOEX/LABEL  UTILIZATION  DEVIATION  UTIL  UTIL  COUNT 

1  SERVICE  TIME  16.4863  8.6246  27  0  2558 


••RESOURCE  STATISTICS** 


RESOURCE 

NUMBER 

RESOURCE 

LABEL 

CURRENT 

CAPACITY 

AVERAGE 

UTIL 

STANOARO  MAXIMUM  CURRENT 
DEVIATION  UTIL  UTIL 

1 

CHECKERS 

1 

16.49 

8.625  27 

0 

RESOURCE 

NUMBER 

RESOURCE 

LABEL 

CURRENT 

AVAILABLE 

AVERAGE  MINIMUM 

AVAILABLE  AVAILABLE 

MAXIMUM 

AVAILABLE 

1 

CHECKERS 

1 

0.1887  -11 

5 

Case  2;  11  Jun  87 

SLAM  II  SUMMARY  REPORT 


SinULRTIOfi  PROJECT  Con  Laada  2  BV  CRPT  FREV 

DATE  11/16/1988  RUM  NUMBER  50  OF  50 


CURRENT  TIME  0.9000E+03 

STATISTICAL  ARRAVS  CLEARED  AT  TIME  O.OOOOE+OO 


♦♦STATISTICS  FOR  URRIRBLES  BASED  ON  OBSERVATION** 

MEAN  STANORRO  COEFF.  OF  MINIMUM  MAXIMUM  NO. OF 

VALUE  DEVIATION  VARIATION  VALUE  VALUE  06S 

TIME  IN  SYSTEM  0.121E+02  0.786E+01  0.649E+00  0.159E+01  0. 154E+03  **** 

TIME  IN  QUEUE  0.696E+01  0.759E+01  0. 109E+01  O.OOOE+OO  0. 150E+03  **** 


♦•FILE  STATISTICS** 


FILE  AVERAGE  STANORRO  MAXIMUM  CURRENT  AVERAGE 

NUMBER  LABEL /TYPE  LENGTH  DEVIATION  LENGTH  LENGTH  WAIT  TIME 


1 

2 


AWAIT  5.856  5.943  27 

CALENDAR  18.051  8.452  30 


0  1.858 

1  0  951 


♦♦REGULAR  ACTIVITY  STATISTICS** 

ACTIVITY  AVERAGE  STANDARD  MRXIMUM  CURRENT  ENTITY 

INDEX/LABEL  UTILIZATION  DEVIATION  UTIL  UTIL  COUNT 

1  SERVICE  TIME  16.3874  8.4792  27  0  2522 


•♦RESOURCE  STATISTICS** 

RESOURCE  RESOURCE  CURRENT  AVERAGE  STRNDARO  MAXIMUM  CURRENT 

NUMBER  LABEL  CAPACITY  UTIL  DEVIATION  UTIL  UTIL 

1  CHECKERS  1  16.39  8.479  27  0 


RESOURCE  RESOURCE  CURRENT  AVERAGE  MINIMUM  MRXIMUM 

NUMBER  LABEL  AVAILABLE  AVAILABLE  AVAILABLE  AVAILABLE 


1  CHECKERS 


0.4144  -11  16 
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Case  3;  11  Jun  87 

SLAtl  II  SUMMARY  REPORT 

SIMULATION  PROJECT  Con  Laada  2  BV  CRPT  FREV 

DATE  11/16/1968  RUM  NUMBER  SO  OF  50 


CURRENT  TINE  0.8000E+03 

STATISTICAL  ARRAYS  CLEARED  AT  TINE  O.OOOOE+OO 


STATISTICS  FOR  URRIA6LES  BASED  ON  OBSERURT I  ON** 

MEAN  STANDARO  COEFF.  OF  NININUN  NAXINUN  NO. OF 

UALUE  DEU I  AT  I  ON  UARIATION  UALUE  UALUE  OBS 

TINE  IN  SYSTEM  0. 139E+02  0. 101E+02  0.726E+00  0. 156E+01  0. 162E+03  **** 

TINE  IN  QUEUE  0.872E+01  0.986E+01  0.113E+01  O.OOOE+OO  0. 157E+03  **** 


**FILE  STATISTICS** 

FILE  AUERAGE  STANDARD  NAXINUN  CURRENT  AUERAGE 

NUMBER  LABEL /TYPE  LENGTH  OEUIATION  LENGTH  LENGTH  WAIT  TINE 


1 

AUAIT 

16.298 

14.204 

63 

0 

5.  101 

2 

CALENDAR 

18.463 

8.482 

30 

1 

0.960 

♦♦REGULAR  RCTIUITY  STATISTICS** 

RCTIUITY 

AUERAGE 

STANDARO 

NAXINUN 

CURRENT 

ENTITY 

INOEX /LABEL 

UTILIZATION 

OEUIATION 

UTIL 

UTIL 

COUNT 

1  SERUICE  TINE 

16.4270 

8.4610 

27 

0 

2556 

**RESOURCE  STATISTICS** 


RESOURCE 

NUMBER 

RESOURCE 

LABEL 

CURRENT 

CAPACITY 

AUERAGE  STANDARO  MAXIMUM  CURRENT 
UTIL  OEUIATION  UTIL  UTIL 

1 

CHECKERS 

1 

16.43  8.461  27 

0 

RESOURCE 

NUMBER 

RESOURCE 

LABEL 

CURRENT 

AUA1LABLE 

AUERAGE  NININUN 

AUAILABLE  AUAILABLE 

NAXINUN 

AUAILABLE 

1 

CHECKERS 

1 

0.3357  -11 

15 
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Case  4:  1 1  Jun  87 

SLAM  II  SUMMARY  REPORT 

SIMULATION  PROJECT  Con  Lamda  2  BY  CRPT  FREY 

DflTE  11/16/1988  RUN  NUMBER  50  OF  50 


CURRENT  TIME  0.8000E+03 

STRTISTICRL  ARRAYS  CLEARED  AT  TIME  O.OOOOE+OO 


♦■"STATISTICS  FOR  URRIABLES  BASED  ON  08SERUAT I  ON*" 

MEAN  STANDARD  COEFF.  OF  MINIMUM  MAXIMUM  NO. OF 

UALUE  DEU I AT  I ON  UARIATION  UALUE  UALUE  08S 

TINE  IN  SYSTEM  0. 163E+02  0. 137E+02  0.838E+00  0. 156E+01  0. 170E+03  **** 

TIME  IN  QUEUE  0.  1 1 1E+02  0.  135E+02  0. 121E+01  O.OOOE+OO  0.  166E+03  **** 


•♦FILE  STATISTICS** 


FILE 

NUMBER  LABEL/TYPE 


AUERAOE  STANORRO  MAXIMUM  CURRENT  AUERAGE 

LENGTH  DEU I  AT  I ON  LENGTH  LENGTH  WAIT  TIME 


1 

2 


AWAIT  10.989  12.032  54 

CALENDAR  17.826  8.451  30 


3.579 

0.964 


♦♦REGULAR  RCTIUITY  STATISTICS** 

ACTIUITY  AUERAGE  STANORRO  MAXIMUM  CURRENT  ENTITY 

INDEX/LABEL  UTILIZATION  DEUIATION  UTIL  UTIL  COUNT 

1  SERUICE  TIME  16.2872  8.4402  27  0  2456 


•♦RESOURCE  STATISTICS** 


RESOURCE 

NUMBER 

RESOURCE 

LABEL 

CURRENT 

CAPACITY 

AUERAGE  STANDARO  MAXIMUM  CURRENT 
UTIL  DEUIATION  UTIL  UTIL 

1 

CHECKERS 

1 

16.29 

8.440  27 

0 

RESOURCE 

NUMBER 

RESOURCE 

LABEL 

CURRENT 

AUAILABLE 

AUERAGE 

AUAILABLE 

MINIMUM 

AUAILABLE 

MAXIMUM 

AUAILABLE 

1 

CHECKERS 

1 

0.4695 

-11 

17 
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